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ABSTRACT 

Infrared (IR) emission spectra are calculated for dust heated by starlight, for mix- 
tures of amorphous silicate and graphitic grains, including varying amounts of polycyclic 
aromatic hydrocarbon (PAH) particles. The models are constrained to reproduce the 
average Milky Way extinction curve. 

The calculations include the effects of single-photon heating. Updated IR absorp- 
tion properties for the PAHs are presented, that are consistent with observed emission 
spectra, including those newly obtained by Spitzer Space Telescope. We find a size dis- 
tribution for the PAHs giving emission band ratios consistent with the observed spectra 
of the Milky Way and other galaxies. 

Emission spectra are presented for a wide range of starlight intensities. We calculate 
how the efficiency of emission into different IR bands depends on PAH size; the strong 
7.7 /xm emission feature is produced mainly by PAH particles containing < 10 3 C atoms. 
We also calculate how the emission spectrum depends on U, the starlight intensity 
relative to the local interstellar radiation field. The submm and far-infrared emission is 
compared to the observed emission from the local interstellar medium. 

Using a simple distribution function, we calculate the emission spectrum for dust 
heated by a distribution of starlight intensities, such as occurs within galaxies. The 
models are parameterized by the PAH mass fraction (/pah, the lower cutoff U m \ n of the 
starlight intensity distribution, and the fraction 7 of the dust heated by starlight with 
U > U m i n . We present graphical procedures using Spitzer IRAC and MIPS photometry 
to estimate the parameters qpAH, U m { n , and 7, the fraction /pdr of the dust luminosity 
coming from photodissociation regions with U > 100, and the total dust mass Md us t- 
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Subject headings: radiation mechanisms: thermal; ISM: dust, extinction; infrared: 
galaxies; infrared: ISM 



1. Introduction 



Interstellar dust in galaxies absorbs energy from starlight; this absorbed energy is then rera- 
diated at i nfrared (IR) and fa r-IR wavelengths. The unprecedented sensitivity of Spitzer Space 
Telescope (jWerner et al.ll2004al ) allows this IR emission to be measured for a wide range of galaxy 
sizes and morphologies. The spectral properties of the IR emission from dust allow one to infer 
the composition of the dust, the size distribution of the dust particles, the intensity of the starlight 
that is heating the dust, and the total mass of dust. 

Deducing the dust properties and starlight intensity is by no means direct or straightforward. 
To start, one needs to assume a provisional dust model with the physical nature of the dust 
(composition, geometry, and size distribution) fully specified. One then tries to constrain the 
grain properties by comparing model predictions with observations (e.g. interstellar extinction, 
scattering, polarization, IR and microwave emission, and interstellar depletions). 

Various models have been propo sed for interstellar dust. The models fall into three b road 



categories: the silicate-graphite model (jMathis et al.lll977l:lDraine fc Lee 



its natural extension - the silicate-gra phite- PAH model (jSiebenmorgen fc Kriigellll992l ; 



2001a 



1984 



Kim et al.lll994 N ) and 



Weingar tner &: Draind l2001aT): the sili cate core-carbonaceous mantle model (jDesert et al 



Li&Draine 



1990j; I Jones et al.lll990l ; iLi &; Greenberglll997l ). and the composite model which assumes the dust 



to be low-density aggregates of small silicates and carbonaceous particles (jMathis fc Whiffen 



Mathis 



1996 



1989 



Zubko et al.l 12004) . The core- mantle model is challenged by the nondetection of the 



3.4 /xm C-H aliphatic hydro carbon polarization on sightlines wh ere the 9. 7 urn Si-Q silica t e ban d 



is observed to be pola rized (jChiar et al 



2006 



Mason et al.1 12006) [se e also 



Adamson et al 



dl999|) 



Li fc Greenbergl (J2002J)]. The original composite model proposed by lMathis Whiffed (|1989l ). with 
~ 80% of the grain volume consist ing of vacuum, may have too flat a far-IR emissivit y to explain the 
observational data (|Draindll994l ), The updated version of the composite model of iMathid (|1996l ) 
assumes a vacuum fraction ~ 45% in order to make economical use of the available carbon to account 
for the observed e xtinc tion while s atisfying the "subsolar" interstellar abundance budget suggested 
bv lSnow k, Wittl (|1996l ). However, iDwek! (|1997l ) argues that this dust model emits too much in the 
far-IR in comparis on with observations. The latest version of the composite model put forward by 
Zubko et al.l (120041 ) aims at reproducing the observed interstellar extinction and IR emission within 
a "subsolar" abundance budget. The ability of composite grain models to reprodu ce the observed 
extinction within "subsolar" abundance constraints was recently challenged by|LiJ (120051 ) based on 
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the Kramers-Kronig relations^ 



In this paper we calculate the IR emission expected from a specific physical model of interstellar 
dust - the silicate-graphite-PAH model. In this model, the dust is assumed to consist of a mixture 
of carbonaceous grains and amorphous silicate grains, with size distributi ons that are consistent 
with t he observed wavelength-dependent extinction in the local Milky Way (jWeingartner fc Draine 
2001al ). including different amounts of polycyclic aromatic hydrocarbon (PAH) material. 



Sp ectroscopic ob s ervat io ns with ISO reveale d the rich spectrum of PAHs (see, e.g.. lBeintema et al 

jlQQfih : iTielens et all jl99fll ): l-Toblin et al.l ||200gI )). The model presented here posits "astro-PAH" 
absorption properties th at are consistent w ith spectroscopic observations of PAH emission from 
dust in nearby galaxies (ISmith et al.l 120061 ) . The model calculations show how the IR emission 
depends on the PAH abundance, and also on the intensity of starlight heating the dust. Compari- 
son of models to infrared photometry obtained with Spitzer Space Telescope allows one to estimate 
PAH abundances, starlight intensities, and total dust masses. 



by 



The plan of this paper is as follows: 

With new labo r atory d ata [particularly the near-IR absorption spectra of PAH ions measured 



Mattioda et al 



( Smith et al 



2004a 



(200 5bf)1 and Spitzer spectroscopic dat a [with new PAH features discovered 



Werner et al. 2004b: Smith et al 



2006) and the familiar 7.7, 8.6, and 11.3 /jm 



features resolved into several subfeatures (jSmith et al.l 120061) we update in ^2] t he cross sections 
C a bs (A) that we previously assumed for the PAH particleJ^I ( Li &: Draine 2001a ). We distinguish 
the PAHs by 2 charge states: either neutral (PAH°) or ionized (PAH 1 * 1 ). We assume that multiply- 
charged (both positively and negatively) PAHs all have the same cross sections as those of singly- 
charged cations. The resulting "astro-PAH" absorption cross sections do not represent any specific 
material, but should approximate the actual absorption properties of the PAH mixture in interstellar 
space. 

In £j3] we show examples of temperature distribution functions for both neutral and ionized 
PAHs heated by starlight. These temperature distribution functions are then used to calculate 
emission spectra for individual particles in £j4j which are in turn used to prepare a plot (Figure [6]) 
showing how the efficiency of emitting into different emission features depends on the PAH size. 



1 Total (gas + dust) interstellar a bundances remain un certain. The most recent determination of the solar oxygen 
abundance (O/H)© = iq~ 3 - 34±0 05 |Asplund et al. 2004 ) is close to the value for B stars in the O rion association 
O/H = io- 3 - 31±0 03 (jCunhaet alJbood ). The warm interstellar medium gas has O/H = lO" 3 - 41 * 01 jCartledge et al ' 
2006), consistent with ~20% of the oxygen in silicate grains. For "refractory" elements such as Mg, Si, and Fe, how- 



ever, partial separation of dust and gas by processes such as radiation pressure, ambipolar diffusion, and gravitational 
sedimentation in the proce ss of star formation could result in stellar abundances differing from interstellar abundances 
bnow 2000; iDraind 12004 ). Photospheric abundances may therefore not be a reliable indicator of interstellar abun- 
dances for the elements that form refractory grains. 



Note that "PAH particles" , "PAH grains" , and "PAH molecules" are synonymous. 
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When we observe emission from a region in the Milky Way or another galaxy, there will always 
be a mixture of dust types and sizes present. We consider the emission from a specific set of 
dust mixtures, as described in £}5j All the dust mixtures considered here are consistent with the 
observed "average" extinction curve for diffuse regions in the local Milky Way, but they differ from 
one another in the assumed abundance of small PAH particles. In £}6] we show that the far-IR 
and submm emission calculated for the model is consistent with the COBE-FIRAS observations of 
emission from dust in the local Milky Way. In $7] we show that the calculated IRAC band ratios 
ap pear to be co n sisten t with the spectrum of diffuse emission from the Milky Way, as determined 



bv lFlagev et al.l (|2006l ) from Spitzer Space Telescope observations. 



The long wavelength emission from the dust model depends on the intensity of the starlight 
heating the dust. It will often be the case that the region observed (e.g., an entire star-forming 
galaxy) will include dust heated by a wide range of starlight intensities. In $8] we describe a simple 
parametric model for the distribution of the dust mass between regions with starlight intensities 
ranging from low to very high. We show in §[9] how observations in the 3 MIPS bands (24 /im, 71 /mi, 
and 160 /im) plus the 3.6 fim and 7.9 /mi IRAC banda^l can be used to estimate the parameters 
describing the distribution of starlight intensities as well as the fraction (/pah of the total dust mass 
that is in PAHs, thereby allowing estimation of the total dust mass Md uat in the emitting region. 
A reader interested primarily in applying the results of this paper to interpretation of IRAC and 
MIPS observations may wish to proceed directly to 53 

The results of the paper are discussed in flOl and summarized in §111 



2. PAH Cross Sections: Post-Spitzer Era 

The strong and ubiquitous interstellar emission features observed at 3.3 /im, 6.2 /mi, 7.7 /mi, 
8.6 /im, and 11.3 /im almost certainly arise from vibrational modes of polycyclic aromatic hydrocar- 
bon (PAH) material, with C-H stretching modes producing the 3.3 /jm feature, C-H bending modes 
producing the in-plane 8.6 /zm and out-of-plane 11.3 /im features, and C-C st retching and bending 



modes producing emissio n features at 6.2 /im, 7.7 /tm, and longer wavelengths (ILeger fc Pugetjll98J ; 
Allamandola et al.lll98.^ FI 



3 The IRAC and MIPS bands are often referred to using nominal wavelengths 3.6, 4.5, 5.8, 8.0, 24, 70, and 160 /tm. 
However, to two significant digits, the effective wavelengths of IRAC bands 3 and 4 are 5.7 and 7.9 /xm (IRAC Data 
Handbook Version 3.0), and that of MIPS band 2 is 71 /im (MIPS Data Handbook Version 3.2) and we will use these 
wavelengths to refer to the bands. 



4 Other carriers have also been propose d, including hydrog enated amorphous carbon ( HAC) ( Jones et al. 1990T I. 
quenched carbo naceous composite CO CO dSakata et al.l 1987J). co al ( Papoufar et alJll993h . hydrogenated fullerenes 
(Webster 1993), and nanodiamonds ( Jones fc d'Hendecourtl 200ol ). The HAC, QCC, and coal hypotheses assume 



that the emission arises following photon absorpt ion in small thermally-isolated aromatic units within or attached 
to these bulk materials (jDulev fc Williams! Il981f ). However, it does not appear possible to confine the absorbed 
photon energy within these aromatic "islands" for the time > 10 -3 s required for the thermal emission process (see 



5 




X (/i.m) A (yu.m) 



Fig. 1. — Observed 5 -20/mi spectra f or: (a) Reflection ne bula NGC 7023 (|Werner et alJbotMbl); (b) Orion Bar PDR 



dVerstraete et al.ll200ll ); (c) M17 PDR JPeeters et al.ll2005h ; (d) Planetary nebula NGC 7027 



20041 ); (e) Sevfert Galaxy NGC 5194 (jSmith et alj|2006l ' 



van Diedcnhoven ct al 



dust model with Qpah 



Also shown (f) is the emission calculated for the present 
4.6%, illuminated by the local diffuse starlight with (7 = 1 and 10 5 (see Fig. I13p . 



To reproduce this emission, a dust model must include a substantial population of ultrasmall 
grains or large molecules with the vibrational properties of PAH material and with sizes such that 
single-photon heating can excite the observed vibrational emission. Because the exact composition 
of the interstellar PAH material is unknown, and also because laboratory knowledge of the optical 
properties of PAHs is very limited, it is necessary to make assumptions regarding the absorption 
cross sections of the PAH particles. 

The approach taken here to modeling the PAHs is to try to find optical properties for "astro- 



Li fc Drainel (2002a)). Free-flying fullerenes and nanodiam onds have the requir ed small heat capacity, but: (1) There 



are strong upper limits on the abundance of C60 and Cg~ (|Moutou et al.lll999h : (2) Although not ruled out, there 
little spectroscopic evidence for hydrogenated nanodiamond in the ISM. 



6 



Table 1: PAH Resonance Parameters. 

CTintj = / CT a bsjrfA 
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PAH" material that appear to be physically reasonable (i.e., band strengths within the range 
measured for PAHs in the laboratory), in order to estimate the sizes and abundance of interstellar 
PAHs that would be consistent with the observed emission spectra. To this end, we adopt feature 
profiles that are based on astronomical observations. A t wavelengths A > the spectra for the 

central regions of galaxies in the SINGS galaxy sample (jSmith et al.ll2006l ) pro vide observed profile s 
for the spectral features that the present model attempts to mimic. We follow ILi Draind (|2001al ) 
(herafter LD01) in describing the PAH vibrational resonances by Drude profiles (components j = 6- 
26 in Tabled])! 



' Boulanger et al. ( 1998h have shown that Lorentz profiles provide good fits to the PAH emission features from the 
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25 



AC abs (A) = A C ]T 



(i) 



but with new profile parameters, adjusted to closely resemble the observed profiles. The profile 
parameters are given in Table dj where Xj, "fjXj, and a- mt = / cr a b s dA _1 are, respectively, the 
peak wavelength, the FWHM, and the integrated strength per C atom of the j-th Drude reso- 
nance. Thefeatures in Table ED and the ir identifications, have been discussed previously [see, e.g., 
Tielend ( 20051 )]. with lSmith et al.l (|2006l ) providing a comprehensive study of the 5-35/mi emission 
from galaxies. We assume the PAH p articles to have ab sorption cross sections per carbon atom 
Cabs(A)/iVc equal to those adopted bv lLi Draind (|2001al ) [see eq. 5-11 of that paper] except for 
the following changes: 



1. For PA H ions, we add additional absorption in the near-IR as recommended bv lMattioda et al 
(|2005al ). This consists of a "continuum" term 

AC abs (A) 



N c 



3.5 x kt 19 - 1 - 45 /* exp (-O.lx 2 ) cm 2 for ions, x = (A//tm) 



-1 



(2) 



Because PAH ions were already assumed to absorb strongly at A < 0.8/mi, this additional 
absorption is numerically insignificant for A < 0.8/mi; the factor exp (— O.lx 2 ) has been added 
simply to force this term to go smoothly to zero as A — > 0. 

2. For PAH ions, we add near-IR resonances at wavelengths 1.05/mi and 1.26/mi plus a negative 
"r esonance" term at L905 um to suppress aborption in the 1.8-2.0 /mi region, as recommended 
bv lMattioda et al.l (|2005al ); the three features are represented by Drude profiles, with param- 
eters as given in Table [p The negative term at A = 1.905 /mi (J = 5) was suggested by the 
recent laboratory data of iMattioda et al.l (|2005bl ). Removal of this term has a negligible ef- 
fect on the heating or cooling rates of PAHs except in regions illuminated by very cool stars 
(T e g < 1500 K) or in regions where the PAHs are excited to unusually high temperatures 
(T> 1500 K). 

3. Small changes have been made to central wavelengths and feature widths (e.g., Xj = 6.20/iin — ► 
6.22/mi and jj = 0.032 — > 0.0284 for the feature near 6.2 /mi, Xj = 11.9/mi — ► 11.99/xm 
and 7j = .025 — > 0.050 for the feature near 12 /mi, and A,- = 12.7/zm — ► 12.61/im and 
7j = .024 — > 0.0435 for the feature near 12. 7 /mi) guided by spectra obtained recently by 
Spitzer Space Telescope (jSmith et al.l 120061 ) . 



4. The integrated strength a- m t = j c^dX of the 3.3 /mi feature has been increased by a factor 
1.5 for neutrals, and a factor 2 for ions, to better agree with the range of values calculated 
for a number of PAHs (see Figure [2|) . 



NGC 7023 reflection nebula and the p Oph molecular cloud. The Drude profile closely resembles a Lorentz profile 
(both have more extended wings than a Gaussian profile). We favor the Drude profile as it is expected for classical 
damped harmonic oscillators. 



s 



5. o"i n t for the 6.22/im feature is 50% of the value in LD01. 



6. The 7.7 /im complex is now composed of three components, at 7.417, 7.598, and 7.850 /xm, 
with o"i n t equal to 50% of the 7.7 pm feature in LD01. 

7. dint of the 8.6 /im feature in LD01 is now shared by features at 8.330 and 8.610 /im, with a m t 
equal to 50% of the 8.6/im feature in LD01. 

8. The 11.3 fj,m feature is now composed of features at 11.23 and 11.30 /mi, with a mt equal to 
50% of the 11.3 pm feature in LD01. 

9. The integrated strength <7i n t of the 12.7/im feature has been multiplied by 0.63 for both 
neutrals and ions relative to LD01. 

10. Weak features have been added at 5.70 /xm, 6.69 /im, 13.60 /im, 14.19 / im, 15.90 ^m, an d 



18.92 /im, as seen in spectra of star-forming galaxies in the SINGS survey (ISmith et al 



The 5.70/m i feature has previously been seen in planet ary nebulae (jAllamandola et al 



and PDRs (Verstraete et al 



1996 



2006). 



1989al ) 



Peeters et al.ll2004l ). and is presumed to be due to combi- 



nation and overtone bands involving C-H out-of-plane bending modes. 



11. A we a k feature at 5.25/im seen in spectra of the M17 PDR and the Orion Bar (jVerstraete et al 



19961 : iPeeters et al.l 120041) and presumed to be due to C-H out-of-plane combination and 
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overtone modes (jAllamandola et al.lll989al ). has been added. 

The strength of the 16.4 /im feature has been multiplied by 0.14 relative to LD01. 



13. A new emission complex near 17 pm has been added, composed of features at 17.038, 17.377, 
and 17.873 /im dSmith et al.ll2004al : Iwerner et al.l \20MU : I Smith et al.lbood ). 



14. Emis sion features at 21.2 and 23.1 /im were seen in some laboratory samples (jMoutou et al 



19961 ). and were therefore included by LD01 as examples of fea tures that might be observed 



15. 



at A > 20/iin. However, the SINGS spectra ([Smith et al.ll2006l ) do not show any features at 
A > 19/im. The 21.2 and 23.1 /im features have therefore been eliminated in the new model. 

LD01 included a broad absorption component with Xj = 26/im, jj = 0.69, and o"i n t = 
18 x 10~ 20 cm/C. This has been replaced by a broad absorption component with A27 = 15/im, 
727 = 0.8, and a- mt = 50 x 10 -20 cm/C to provide continuum emission from 13/im longward. 



In general, PAHs and larger grains will not be spherical, but we will characterize a grain of mass 
M by the effective radius a, defined to be the radius of an equal volume sphere: a = (3M/47T/)) 1 / 3 , 
where amorphous silicate is assumed to have a mass density p = 3.5 gem -3 , and carbonaceous 
grains are assumed to have a mass density due to graphitic carbon alone of p = 2.2 gem -3 . 



Thus the number of carbon atoms in a carbonaceous grain is 

* 3 

Nq = 460 ' 



10 A 



(3) 
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The smallest PAH considered in this paper ha s Ng = 20 C atoms (corr esponding to a ~ 3.55 A) 
since smaller PAHs are photolytically unstable (|Allamandola et al.lll98 9b). As in LD01, we assume 
the number of H atoms per C atom to depend on the size of the PAH: 



H/C = 0.5 

= 0.5(25/A c ) 1/2 
= 0.25 



for N c < 25 

25 < Ac < 100 
A c > 100 



(4) 



o 

\ 
1/2 . 

CO 

U 
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- terrylene .CjoH,, 




7, 8 9 10 



A(yum) 



15 



20 



Fig. 2. — Solid curve: adopted absorption cross section per C from eq. ((5]) with C/H~3.2 (e.g., C64H2o). For the 
neutrals, anions, and cations listed in the figure legends, the horizontal li ne segments i ndicat e the average absorption 
over that frequency intervals, taken from theoretical calculations by Malloci et al. 1 20061 ). The heavy solid line 
segment is the average for the 6 species shown. 



Figure [2] shows the adopted C a b s (A) in the infrared for a PAH molecule with H/C=5/16 (e.g., 
C64H20, with a = 5.18 A). Also shown are values of C a b s per C atom, averaged ov er wavelength 
intervals, for a number of molecules for which C a b s has been calculated theoretically (jMalloci et al. 
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2006) for selected PAH molecules, cations, and anions. The first thing to note is the wide range of 
absorption cross sections per C. For example, in the case of the 3.3 fim C-H stretch, the integrated 
absorption cross section per C varies by a factor of 25 among the PAH cations. The absorption 
averaged over the 7.5-8.0 /im range varies by a factor of 30 among the neutral PAHs. Similar large 
variations in absorption cross sections are also seen at other wavelengths. Our adopted cross section 
falls within the range found for the sample of molecules shown in Fig. [2j 

Figure [3] compares our adopted absorption cross sections for both neutral and ionized PAHs. 

As in LD01, as the number Nq of carbon atoms in the grain increases, we assume a continuous 
change in optical properties from those of PAH material when Nq is small, to those of graphite 
when Nq is large. The transition from PAH to graphite is entirely ad-hoc: we take 



Cabs (A) = (1 - Cgra) C abs (PAH, N c ) + £ g ra Cabs (graphite, a) 

where we take the graphite "weight" £ gra to be 

£ gra = 0.01 for a < 50 A (A^ < 5.75 x 10 4 ) 

'50A^ 31 



0.01 + 0.99 



1 



for a > 50 A (i.e. JV C > 5.75 x 10 4 



(5) 

(6) 
(7) 



The rationale for this is as follows: in addition to the C-H stretching mo de emission at 3.3 /xm, there 



appe ars to be 2-5 /im continuum emission from the interstellar medium (|Lu et al 



2003 



Helou et al 



2004) and we therefore need a source of continuum opacity in the 2-5 /xm region that is not provided 
by the C-H and C-C stretching and bending modes. Here we assume that every small PAH has 
a small amount of "continuum" opacity, equal to 1% of what would have been calculated with 
the optical properties of bulk graphite (i.e. £ gra = 0.01). When the carbonaceous particles have 
Nc ^ 10 5 , they are assumed to behave like bulk graphite. 

Graphite is highly anisotropic, with different dielectric functions e± and en for electric fields 
perpendicular and parallel to the "c-axis" (the c-axis is normal to the "basal plane" of graphite). 
For l^m < A < 20/im, absorption in small randomly-oriented graphite spheres is primarily due to 
the free electrons in graphite moving in the basal plane. However, the basal plane conductivity 
is large enough that at wavelengths A > 30/xm the absorption is primarily due to the weak but 
nonzero conductivity parallel to the c-axis. The c ontribution of "free electrons" to en results in an 
absorption peak near 30 jJLm. ( Draine Led Il984j ). This peak is seen for a > 60 A in Figure [3b- 
This absorption peak results in a broad emission feature near 30 //m when the a > 60 A graphite 
particles are heated to T > 100 K. 

The peak in C a b s near 30/mi in Figure [3b is a consequence of our adopted dielectric function for 
graphite, which is based on a simple "free electron" model that may not apply to real istic carbona- 



ceous grains. For example , the measured absorption in amorphous carbon grains (jTanabe et al 



19831 : iMennella et al.lll999l ) does not appear to show a peak near 30/im. This is further discussed 
in TOTTl 
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Fig. 3. — Absorption cross section per C atom for (a) neutral and ionized PAHs, and (b) for ionized carbonaceous 
grains, with properties of PAHs for 6 A < a < 50 A, and properties of graphite spheres for a > 100 A. See ^3 for 
details. 



3. Heating of Dust by Starlight 



We consider heating of grains by radiation with energy density per unit frequency 



U x u 



MMP83 



(8) 



wh ere U is a dimension less scaling factor and u 



MMP83 



is interstellar radiation field (ISRF) estimated 



by iMathis et al.l (119831 ) [hereafter MMP83] for the solar neighborhood. 



For each grai n composition and radius a, we use a detailed model for the heat capacity 
(jDraine Lill200ll ) to calculate the function E(T), the expectation value for the vibrational en- 
ergy of the grain when in equilibrium with a heat reservoir at temperature T. For a grain with 
vibrational energy E, the grain "temperature" T(E) is taken to be the temperature at which the 
expectation value for the vibrational energy would be E: E(T) = E; we use t his temperature est i- 
mate for all values of the vibrational energy E, in this respect departing from lDraine Lil (|200ll ). 
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who used a different estimate for the temperature when dealing with the first 20 vibrational excited 
states. This does not appreciably affect the emission spectrum since almost all of the absorbed 
photon energy is reradiated while the grain is at high temperatures. 

For each grain composition, radius a, and radiation intensity scale factor U, we determine 
the probability distribution function dP/dT, where dP is the probability of finding the grain 
with temperature in [T, T+ dT\. For large grains, dP/dT is approximated by a delta function 
dP/dT 5(T — T ss ), where the "steady-state" temperature T ss (a) is the temperature at which the 
radiated power is equal to the time-averaged heating rate for a grain of radius a. For small grains, 
we find the steady-state solution dP/dT for grains subject to stochastic heating by photon ab- 
sorption, and coo ling by emission of infrared photons, using the "thermal-discrete" approximation 



(jDraine Sz Lill200ll ) , where the downward transition probabilities for a grain with vibrational energy 
E are estimated using a thermal approximation. For each grain size a and radiation intensity U, 
we divide the energy range [E m i n , -E max ] into 500 bins. -Emm and E ma _ x are found iteratively, with 
the requirement that the probability of the grain being outside the range [E m i n , E max \ be negligible. 

Figured] shows dP/dlnT for PAH + /graphite grains for selected grain sizes, for U = 1 and U = 
10 4 . In FigureHK one sees that small grains undergo extreme temperature excursions (the o = 3.55 A 
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PAH occasionally reaches T > 2000 K), whereas larger grains (e.g., a = 300 A) have temperature 
distribution functions that are very strongly-peaked, corresponding to only small excursions around 
a "steady-state" temperature T ss (the temperature for which the rate of radiative cooling would 
equal the time-averaged rate of energy absorption). 

Figure shows dP/dlnT for U = 10 4 . It is apparent that when the rate of starlight heating 
is increased, the "steady-state" temperature approximation becomes valid for smaller grains. For 
example, for U = 10 one could approximate an a = 50 A grain as having a steady temperature 
T ss ~ 150 K whereas for U = 1 the temperature excursions are very important for this grain. The 
radius below which single photon heating is important, and above which the grain temperature 
can be approximated as being constant, is the size for which the time between photon absorptions 
becomes equal to the radiative cooling time (jDraine fc Lill200ll ). Or (equivalently) , it is the size 
for which the thermal energy content of the grain when at T ss is equal to the energy of the most 
energetic photons heating the grain. 



4. Single-Grain Emission Spectra 

From the probability distributions dP/dT, we calculate the time-averaged emission spectra for 
individual particles, 

f dP 
Px = I 4vrC abs (A)B A (r)— dT , (9) 

B X (T) = 2hc 2 X~ 5 [exp(hc/\kT) - . (10) 

Figure [5] shows the 3-30 /im emission from PAH ions and PAH neutrals heated by U = 1 starlight, 
with spectra shown for a number of different sizes. As expected, the short wavelength emission 
(e.g., the 3.3 //m feature) is strong only for the smallest particles, which can be heated to T > 10 3 K 
by absorption of a single ultraviolet photon (see Fig. [4]). As the particle size becomes larger, the 
short wavelength emission falls off, and an increasing fraction of the absorbed starlight energy is 
radiated in the longer wavelength modes. Emission in the 17 ^m complex, for example, is most 
efficient for a ~ 15 A. 

The dependence of the feature emission on grain size is shown in Figure [6l From this plot 
it is apparent that emission in the 3.30 /im feature will be almost entirely due to a < 6 A, or 
Nq < 10 2 , whereas grains with radii as large as 12 A, or Ac ~ 10 3 , are efficient at converting 
absorbed starlight into 7.7 ^m emission. The 11.3 /im feature can be produced by particles as large 
as a ~ 20 A or Nq ~ 4000, and the 16.45 ^m feature and 17.4 fim complex are efficiently produced 
by PAH particles in the 8-25 A size range (300 < N c < 10 4 ). 

The PAH emission spectra shown in Figure [5] and the band emission efficiencies shown in 
Figure [6] are independent of the starlight intensity U in the single-photon-heating regime, U < 10 4 
for A < 30/xm. For larger grains, and at higher starlight intensities, the emission spectra do depend 
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Fig. 5. — Normalized time-averaged emission spectra for U < 10 4 for (a) neutral and (b) ionized PAHs of various 
sizes (see text). 



on U. If C a bs oc A 2 for A > 30/rni, then the power radiated by the grain oc T 6 , the steady-state 
temperature T ss oc U 1 ^, and the infrared emission will peak at a wavelength A p oc 1/T SS oc U^ 1 ^. 

The PAHs in the interstellar medium will consist of a mixture of neutral and ionized particles. 
For a given size PAH, the ionization balanc e will depend on the gas tem perature, the electron 
density, and the ultraviolet rad iation field (Weinga rtner fc Draind l2001bl ). Here we adopt the 
ionization balance estimated by [Li &: Draind (|2001al ) for the diffuse ISM in the Milky Way - a 
weighted sum of the ionization fractions calculated for PAHs in the CNM (cold neutral medium), 
WNM (warm neutral medium) and WIM (warm ionized medium) . For this weighted sum of neutral 
and ionized PAHs, \p\jp at selected wavelengths is shown as a function of grain size in Figure [71 
Note that most of the selected wavelengths coincide with the positions of PAH emission features, 
thus indicating how PAH size affects the efficiency of converting starlight energy into IR emission 
features. The adopted ionization fraction as a function of grain size is shown in Figure [8l 

Figure [7J also shows \p\/p evaluated at the 23.7 p wavelength of the MIPS 24/iin band. We 
see that PAHs in the 15-40A size range (2000 < Nq <3x 10 4 ) are relatively efficient at converting 
starlight energy into 24^m continuum following single-photon heating (Figure H] shows that a single 
hv < 13.6 eV photon can heat a 20 A PAH to ~ 170 K). The term "very small grains" is sometimes 



15 




Fig. 6. — Efficiency for radiating in different emission bands, as function of size, for ionized and neutral PAHs (see 
text). 

used to describe small grains that contribute continuum emission into the IRAS 25/im band or the 
MIPS 24/xm band - here we see that such grains must have effective radii in the 15-40A range to 
efficiently convert absorbed stellar energy into ~24 / um continuum emission. 

Figures [9] and [10] show 2-200 /xm emission spectra calculated for carbonaceous grains and 
amorphous silicate grains with sizes extending from a = 3.55 A to 5000 A. As expected, the 
emission from a > 60 A grains peaks at A 100f7 _1 / 6 / um. However, for U = 10 6 and 10 7 , \p\ for 
a > 60 A carbonaceous grains has two peaks: one near the peak in AB\ near 100f7 _1 / 6 ^m, and the 
second near 30 /xm, the latter resulting from the peak in C a b s (A) near 30 /im (see Figure [3b). As 
discussed in $2] above, the 30/iin opacity peak in the present model produced by graphite in the 
present model may not apply to the carbonaceous material in interstellar dust. 



5. Dust Mixtures 

Interstellar dust in the Milky Way and other galaxies includes a wide range of grain sizes. Here 
we consider the size distributions put forward by WD01 to reproduce the wavelength-dependent 
extinction on Milky Way sightlines with extinction curves characterized by Ry = Ay /{Ab — Ay) = 
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Fig. 7. — Xp\ divided by the the total time-averaged power p radiated by the grain, as a function of grain size 
a, for selected values of wavelength A; curves are labelled by the value of A(/xm). We assume a mix of neutral and 
ionized PAHs, with ionization fraction as in Figure [8] 24/im is the nominal wavelength of the MIPS band; the other 
wavelengths are at peaks of emission features. 



3.1. WD01 included a pop ulation of small (a < 50 A) carbonaceous particles with specified total 



mass. 



Li fc Draind (|2001al ) showed that if these small carbonaceous particles had the physical 



properties of PAHs and were distributed in two log-normal components, 



dn riQj 

T = — ex P 
da ^-^ a 



(\n(a/a Qj )f 
2a) 



+ non — log— normal contribution 



(11) 



the result ing infrared emiss i on was approximat ely consistent with the diffuse emission observed 



by IRTS (lOnaka et al. 



1996 



Tanaka et al.lll996l ). The non-log- normal contribution to eq. (jlip is 



given by eq. (4) of WD01. with parameter values taken from that paper except f or a reduction in 
the numbers of grains per H by a factor 0.92, as recommended by iDraind (|2003l ). Note that the 
non-log-normal term also extends continuously down to the smallest sizes. 

We continue to use eq. (fTTI) for dn/da, but we have modified the values of the parameters aoj 
and <7j, as given in Table EJ The factors riQj in eq. (jlip are related to the numbers bj of carbon 
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Fig. 8. — Adopted fractional ionization in diffuse ISM ( Li fc Draine 2001ah 



Table 2: PAH Size Distribution Parameters 



parameter 


LD01 


present paper 


61/(61 +62) 


0.75 


0.75 


62/(61 + 62) 


0.25 


0.25 


aox(A) 


3.5 


4.0 


ao2( A) 


30. 


20. 


<Ji 


0.4 


0.4 


CT2 


0.4 


0.55 



atoms per total H in each of the log-normal components: 



3 exp(4.5a 2 ) , /HJ 



mc 

( 27r )3/2 1 + erf p^-o 



Mj u 3 

x _ M a Mj/«min) ^ 

where pc is the carbon mass density, mc is the mass of a carbon atom, and 

auj = a 0j exp (3cr|) (14) 

is the location of the peak in the mass distribution oc a?dn/d In a. Note that the non-log-normal 
term in eq. (jlip also contributes to the population of Nq < 10 3 PAHs. 

These size distributions have been constructed for various amounts of carbonaceous material 
in the very small PAH particles. The size distributions all reproduce the observed wavelength- 
dependent extinction in the Milky Way, for sightlines with Ry = Ay/E(B — V) ~ 3.1. Because the 
abundance of ultrasmall grains is important for the infrared emission as well as other applications, 
Table [3] gives the amount of carbon present in grains containing < 10 2 , < 200, < 500, < 10 3 , 



18 




X(fim) X(fim) A(jLtm) 



Fig. 9. — Emission for selected sizes of carbonaceous grains for U = 0.1 - 10 7 . For a < 100 A PAHs are assumed to 
have Xion from Fig. [8] (see text). 



< 10 4 , and < 10 5 C atoms, for the 7 different grains models. Models jm = 1 — 7 have (&i + 62) = 
0.92 x 60 x 10" 6 x [(j M - l)/6]. The mass distributio ns for carbonaceous grains are shown in Fig. 
[TT1 [The silicate mass distribution is shown in Fig. 2 of lWeingartner fc Draind (|2001al ). except that 
the dust abundances should be multiplied by 0.92.] 



The emissivity per H nucleon for a dust mixture heated by starlight intensity U is 



y [dap- 
^ . da 



C a bs{j,a.,v)B v (T) 



(-) 



dT 



(15) 



V 



Bu(T) 



2hv 3 



1 



(16) 



c 2 exp(his/kT) — 1 ' 

where the sum is over compositions j, and the temperature distribution function dP/dT depends 
on composition j, radius a, and starlight intensity U. 
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Fig. 10. — Emission for amorphous silicate grains of various sizes for [7 = 0.1 - 10 7 (see text). 



Taking dP/dT calculated for the local starlight intensity U = 1, summing over the grain size 
distribution, and over the fractional ionization shown in Fig. [8], we obtain the emission per H nucleon 
j\ shown in Figure [T2l The emission at A < 20/xm depends strongly on the PAH abundance (/pah- 

Figure [13] shows the emission calculated for the model with <?pah = 4.6%, but for different 
starlight intensities U. For U < 10 3 , the normalized emission at A < 20/im is essentially independent 
of U, because the emission is almost exclusively the result of single-photon heating, with the PAH 
particles cooling off almost completely between photon absorptions. This remarkable invariance 
of the PAH emission spectra over orders of mag nitude variations i n star l ight intensities has bee n 
obse rved in a wid e rang e of environments [e.g. see iBoulanger et al.l (|2000l ) . iKahanpaa et al.l (|2003l ). 
and ISakon et al.1 (|200J)]. For U > 10 4 , however, the mean time between photon absorptions 
becomes shorter than the radiative cooling time for a PAH with ~1 eV of internal energy, so that 
the small grains do not cool completely between photon absorptions. Photon absorptions are then 
able to take them to higher peak temperatures, and the fraction of the power radiated at A < 20/um 
increases with increasing U. 
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6. Far-IR and Submm Emission 

The focus in the present dust model has been on the PAH features, but the model is intended 
to reproduce the thermal dust emission at far-infrared (FIR) and submm wavelengths as well. 
The opacities of the amorphous silicate and graphite particles in the model are calculated using 
the dielectric functions. The graphite dielectric function is based on laboratory measurements 
of graphite. Because the nature of interstellar amorphous silicate m aterial is unce r tain, the far 



infrared and submm behavior is not constrained by laboratory data. iDraine &: Led (|1984 ) made 



an estimate of the imaginary part of the dielectric function of interstellar amorphous silicate. 



Finkbeiner et al.1 (ll999l . hereafter FDS) used COBE-FIRAS observations of the sky at high 
galactic latitudes, after removal of the cosmic background radiation and zodiacal emission, to 
characterize the emission from diffuse gas and dust in the Milky Way. They excluded |6| < 7°, 
the Magellanic Clouds, HII regions in Orion and Ophiuchus, and an additional 16.3% of the sky 
where the data were of lower quality. The final data set comprises 81% of the |6| > 7° sky. For 
this region, FDS find an empirical fit that quite accurately reproduces the observed 100^m-3mm 
spectra, using two parameters for each pixel: the 100/im intensity I uo and a temperature T 2 that 
determines the shape of the A > 100/xm spectrum: 

\2.70 d (rp \ , n cn e/,,a, \1.67d fm.\ f tl \ L182 



(^ ) 2 - 70 ^(T 2 ) + 0.515(^ ) 1 - 67 ^(r 1 ) / T 2 V 

" Uo B U0 (T 2 ) + 0.515B„(Ti) ' 1 ' 



(17) 
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Fig. 12. — Emission spectra for Milky Way dust models with various gpAH (dust models jj^=\- 
7), heated by starlight with U = 1. Tabulated spectra for this and other cases will be available at 
|http://www.astro.princeton.edu/~draine/dust/irem .html 



where vq = c/100/im = 3000 GHz. The mean value (T2) = 16.2 K. Variations in T2 are presumably 
the result of variations in the intensity of the starlight heating the grains; a simple model suggests 



16.2 KK 



1/6.70 



FDS 



9.4 KU, 



1/5.67 
FDS 



(18) 



where ?7fds is the intensity of the starlight heating the dust, relative to the average for the region 
analyzed by FDS. 



Based on the observed spectrum of eq. (|17p with T2 = 16.2K. lLi &: Draind (|2001al ) made small 
adjustments to the imaginary part of the dielectric function for the amorphous silicate material 
at A > 250/im to improved agreement with the observed emission spectrum. The real part of 
the dielectric functi on is obtained from the imaginary part using the Kramers-Kronig relations 
(jDraine &; Ledll984l ). The resulting dielectric function is used to calculate absorption cross sections 
for amorphous silicate spheres in the present work. 

The upper panel of Figure [FIr shows the observed far- infrared and submm emission spectrum, 
as given by eq. (fTTJ) with T2 = 16.2 K (i.e., t/pDS = 1)- Also shown are emission spectra calculated 
for the present model for two values of the starlight intensity: U = 0.8 and U = 1. Both spectra 
have shapes that are close to the observed spectrum, but the U = 0.8 spectrum agrees to within 
a few percent from 100/mi to 1.5 mm. Evidently C/fds ~ 1.2J7 gives the correspondence between 
the starlight inte nsity U in the presen t model (normalized to the estimate of the local interstellar 
radiation field by iMathis et al.l (|1983l )) and the starlight intensity C/fds relative to the "average" 
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Fig. 13. — Emission spectra for size distribution jm = 7 for selected starlight intensity scale factors U. 

for the 71% of the sky analyzed by FDS. 

Fig. 114b shows the emission spectra for the model when the radiation field is lowered to 
U = 0.5 and raised to XJ = 2; these spectra compare well with the FDS spectra for £/fds = 0.6 and 
[/peg = 2.4, as expected from the correspondence C/fds ~ 1-2C7 inferred from Fig. [T4"a, 

We conclude that the present dust model can successfully reproduce the emission observed 
from dust in the diffuse interstellar medium of the Milky Way (including dust in molecular clouds 
at |6| > 7deg and away from Orion and Ophiuchus) out to wavelengths as long as 2mm without 
introduction of additional emission components. 



23 



Table 3: Physical Dust Models 

C/H (ppm) in PAHs with 



3m 


Model 


<7pah(%) 


(M dust /M H r 


N c < 100 


iV c < 200 


iV c < 500 


N c < 10 3 


iV c < 10 4 


N c < 10 5 


1 


MW3.1.00 


0.47 


0.0100 


1.2 


1.8 


2.9 


3.9 


8.7 


17.0 


2 


MW3.1.10 


1.12 


0.0100 


3.5 


5.6 


8.0 


9.4 


13.6 


21.5 


3 


MW3.1.20 


1.77 


0.0101 


5.8 


9.4 


13.2 


14.9 


18.6 


26.0 


4 


MW3.1.30 


2.50 


0.0102 


8.3 


13.4 


18.6 


20.8 


24.2 


31.1 


5 


MW3.1.40 


3.19 


0.0102 


10.9 


17.6 


24.4 


27.1 


30.7 


37.9 


6 


MW3.1.50 


3.90 


0.0103 


13.5 


21.8 


30.2 


33.4 


37.1 


44.4 


7 


MW3.1.60 


4.58 


0.0104 


16.1 


26.1 


36.1 


39.8 


43.9 


51.6 


8 


LMC2.00 


0.75 


0.00343 


0.9 


1.2 


1.8 


2.2 


3.5 


5.0 


9 


LMC2.05 


1.49 


0.00344 


1.6 


2.6 


3.7 


4.3 


5.2 


6.7 


10 


LMC2.10 


2.37 


0.00359 


2.6 


4.4 


6.3 


7.1 


7.9 


9.6 


11 


SMC 


0.10 


0.00206 


0.1 


0.1 


0.2 


0.3 


0.4 


0.5 



a M dust /M gas = (1/1.36) (M dust /M H ) 



7. Spitzer IRAC and MIPS Band Ratios 

For interpreting observations with the IRAC and MIPS cameras on Spitzer Space Telescope, 
the quantities of interest are the band-convolved emissivities and luminosities 

(j ) b d = J Rba,nd(v) judl> . ^ _ / R hand (v)L v dl> 

JV/^band^-RbandH^ ' " ^ J(^/^band) /3 -Rband(^)^ 



and 



(^Jy)band = ^and ( ju ) band , ( z/ -^y)band = ^and (-^i/)band , (20) 



where i^bandl^) is the relative response per unit power for the combination of optics, filter, and 
detector and i^band = c/Aband) where Aband the nominal wavelength of the band, is given in Table [H 

Table H] gives (i / j'i/)band fo r the four IRAC bands@ the 2 "Peakup" bands of the Infrared Spec- 
trograph (IRS)0 the three MI PS bands§. Table [ 5] gives emissivities convolved with photometric 



bands of the AKARI satellite ( Kawada et al. 2004l ; lonaka et aL 2004 ). and model emissivities con 



volved with the the three bands of the Herschel PACS instrument In, and the three bands of the 



6 iiband(^) oc {l/hv)Shimd {v) , where the relative response per photon S'band(^) is obtained from 
http://ssc.spitzer.caltech.edu/irac/spectral_response.html. The IRAC calibration uses /? = — 1, 

7 -Rband(^) oc ( 1 / hv) 5*band {v) , where the relative quantum efficienty Sband(^) is given in Fig. 6.1 of the Infrared 
Spectrograph Data Handbook Version 2.0. The IRS Peakup calibration uses /3 = — 1. 

^band(^) is obtained from http://ssc.spitzer.caltech.edu/mips/spectral_response. html] The MIPS cali- 
bration is for a 10 4 K blackbody, i.e., /3 = 2. 

9 A. Poglitsch has kindly provided provisional i?band(i / ) for PACS. The calibration procedures for PACS have not 
yet been finalized. Here we assume (3 — 2. 
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Fig. 14. — (a) Squares: observed spectrum, relative to the 100/xm emission, of Milky Way dust from lFinkbeiner et al.l 
(1993), given by eq. (|17I18|I with T2 = 16.2 K. Solid curve: 100//m-3mm emission spectrum for MW dust model with 
q PAH = 4.6%, for U = 0.8 and 1. The U — 0.8 model approximately reproduces the observed emission for T2 = 16.2 K, 
though having excess emission at A > 1.5 mm. The good agreement between the U = 0.8 model and the observed 
spectrum with (7fds = 1 suggests a correspondence C/fds ~ 1.2(7. Dotted curve: single-temperature dust model 
with opacity k cx v 1 and T = 18 K, for comparison. A single modified blackbody produces insufficient emission for 
A > 500/im. (b) Observed spectra for !7fds = 0.6 and 2.4 (squares) and the present dust model for U = 0.5 and 2 
(solid curves). 



Herschel SPIRE instrument l 10 l ax6 given in Table [6} Emissivitios {i^ju )band are given for selected 
grain models with several different values of qpAH, and for different values of the starlight intensity 
scale factor U. 

Figure [15] shows (^L J/ }b an d relative to the total infrared (TIR) dust luminosity 

/"OO 

L T m= / L v dv (21) 

J 

for 5 bands: IRAC 3.6 /tm, IRAC 7.9 /im, MIPS 24 /im, MIPS 71 fim, and MIPS 160 /im. Results are 



10 M. Griffin has kindly provided provisional i?band(^) for SPIRE. The calibration procedures for SPIRE have not 
yet been finalized. Here we assume /3 = 2. 
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Fig. 15. — IRAC 3.6 ^m, IRAC 7.9 pm, MIPS 24 /xm, MIPS 71 /xm, and MIPS 160 /im band strengths for dust 
illuminated by a single starlight intensity U, as a function of U, for dust models with 7 different values of (/pah. 

given for grain models jm = 1-7, with the model results labelled by the value of (/pah- As expected, 
the 3.6 /im emission and 7.9 /im emission per unit total power is entirely the result of single-photon 
heating at low values of U, and therefore (z^i^.g/^TrR is independent of U for U < 10 3 , and 
(^L^)3.6/Ltir is independent of U for U < 10 5 . (vL^^q/Ltir and {uL v )j^/ Xtir are also both 
approximately proportional to (/pah m the low-C/ limit, as expected. For small values of U, the 
24 /im emission is also the result of single-photon heating, and (vL v )2Al '^tir is independent of ?7 
for U < 10. 



Flagey et al.l (120061 ) extracted the diffuse ISM emission from IRAC imaging of various regions 
of the Milky Way obtained by the Galactic First Look Survey (GFLS) and the Galactic Legacy 
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Table 4. Model Emissivities (^jy)band (ergs s sr /H) a for Spitzer Space Telescope 



u 


model 

<?PAH (%) 


IRAC 1 
3.550/xm 


IRAC 2 
4.493/im 


IRAC 3 
5.731/xm 


IRAC 4 
7.872/im 


IRS PU1 
16/im 


IRS PU2 
22./im 


MIPS 1 
23.68/im 


MIPS 2 
71.42^m 


MIPS 3 
155.9/xm 


0.5 


0.47 


5.66e-28 


1.54e-28 


1.31e-27 


4.05e-27 


4.86e-27 


5.91c-27 


5.83e-27 


5.04e-26 


1.25e-25 


0.5 


1.12 


1.55e-27 


4.30e-28 


3.53e-27 


9.17e-27 


6.41e-27 


6.27e-27 


6.05e-27 


4.76e-26 


1.21e-25 


0.5 


1.77 


2.56e-27 


7.09c-28 


5.77c-27 


1.43e-26 


7.97e-27 


6.62e-27 


6.26e-27 


4.54e-26 


1.19e-25 


0.5 


3.19 


4.78o-27 


1.34e-27 


1.07c-26 


2.60o-26 


1.24e-26 


8.87e-27 


8.160-27 


4.11e-26 


1.12o-25 


0.5 


4.58 


7.07e-27 


1.98e-27 


1.58e-26 


3.81o-26 


1.73o-26 


1.14e-26 


1.04e-26 


3.66e-26 


1.05c-25 


1.0 


0.47 


1.13c-27 


3.08e-28 


2.63c-27 


8.10e-27 


9.84e-27 


1.21e-26 


1.20e-26 


1.35e-25 


2.27e-25 


1.0 


1.12 


3.11e-27 


8.59e-28 


7.06c-27 


1.84e-26 


1.29e-26 


1.29e-26 


1.24e-26 


1.28e-25 


2.21c-25 


1.0 


1.77 


5.11e-27 


1.42e-27 


1.15o-26 


2.87o-26 


1.60e-26 


1.35e-26 


1.28e-26 


1.22e-25 


2.18c-25 


1.0 


3.19 


9.55e-27 


2.68e-27 


2.14e-26 


5.20e-26 


2.50e-26 


1.80e-26 


1.66e-26 


l.lle-25 


2.06e-25 


1.0 


4.58 


1.41e-26 


3.96e-27 


3.15o-26 


7.62e-26 


3.47e-26 


2.32e-26 


2.11e-26 


9.89e-26 


1.94c-25 


3.0 


0.47 


3.39e-27 


9.24e-28 


7.90e-27 


2.44e-26 


3.07e-26 


3.95e-26 


3.92e-26 


5.80e-25 


5.19e-25 


3.0 


1.12 


9.32e-27 


2.58c-27 


2.12e-26 


5.51e-26 


3.99e-26 


4.14e-26 


4.03e-26 


5.53e-25 


5.06e-25 


3.0 


1.77 


1.53e-26 


4.25e-27 


3.47e-26 


8.62e-26 


4.91e-26 


4.32e-26 


4.12e-26 


5.35e-25 


5.01e-25 


3.0 


3.19 


2.86e-26 


8.05e-27 


6.41e-26 


1.56e-25 


7.58e-26 


5.66e-26 


5.25e-26 


4.89e-25 


4.77c-25 


3.0 


4.58 


4.24e-26 


1.19c-26 


9.46e-26 


2.29e-25 


1.05e-25 


7.22e-26 


6.60e-26 


4.41e-25 


4.57e-25 


5.0 


0.47 


5.66e-27 


1.54e-27 


1.32e-26 


4.07e-26 


5.28e-26 


6.97e-26 


6.95e-26 


1.08e-24 


7.30e-25 


5.0 


1.12 


1.55e-26 


4.29e-27 


3.53e-26 


9.19o-26 


6.80e-26 


7.27e-26 


7.10e-26 


1.03e-24 


7.13c-25 


5.0 


1.77 


2.56e-26 


7.08e-27 


5.78e-26 


1.44e-25 


8.32e-26 


7.54e-26 


7.23e-26 


1.00e-24 


7.06e-25 


5.0 


3.19 


4. 77c-26 


1.34e-26 


1.07c-25 


2.60e-25 


1.28e-25 


9.76e-26 


9.10c-26 


9.20e-25 


6.76e-25 


5.0 


4.58 


7.06e-26 


1.98e-26 


1.58e-25 


3.81o-25 


1.77c-25 


1.24e-25 


1.14e-25 


8.36c-25 


6.50e-25 


10.0 


0.47 


1.13o-26 


3.08e-27 


2.63c-26 


8.17c-26 


1.12e-25 


1.54e-25 


1.54e-25 


2.37e-24 


1.13c-24 


10.0 


1.12 


3.11e-26 


8.59e-27 


7.07c-26 


1.84e-25 


1.42e-25 


1.59e-25 


1.56e-25 


2.27c-24 


l.lle-24 


10.0 


1.77 


5.11e-26 


1.42e-26 


1.15c-25 


2.880-25 


1.72e-25 


1.63e-25 


1.57e-25 


2.21e-24 


1.10o-24 


10.0 


3.19 


9.55e-26 


2.68e-26 


2.14e-25 


5.21c-25 


2.61e-25 


2.07o-25 


1.94e-25 


2.05e-24 


1.06e-24 


10.0 


4.58 


1.41e-25 


3.96e-26 


3.15c-25 


7.63e-25 


3.59e-25 


2.59e-25 


2.39e-25 


1.88e-24 


1.02e-24 


100.0 


0.47 


1.13c-25 


3.09e-26 


2.67c-25 


8.67e-25 


1.68e-24 


2.80e-24 


2.89e-24 


2.16e-23 


3.72c-24 


100.0 


1.12 


3.11e-25 


8.62e-26 


7.12c-25 


1.89e-24 


1.94e-24 


2.78o-24 


2.83c-24 


2.09e-23 


3.66c-24 


100.0 


1.77 


5.12e-25 


1.42e-25 


1.16c-24 


2.93e-24 


2.20e-24 


2.75e-24 


2.77e-24 


2.06e-23 


3.64e-24 


100.0 


3.19 


9.57e-25 


2.69e-25 


2.15e-24 


5.26e-24 


3.07c-24 


3.09e-24 


3.03e-24 


1.95e-23 


3.55e-24 


100.0 


4.58 


1.41e-24 


3.97e-25 


3.17e-24 


7.70e-24 


4.07e-24 


3.56e-24 


3.41e-24 


1.85e-23 


3.48e-24 



a Additional models, with spectra, are available online at http://www.astro.princeton.edu/~draine/dust/irem.html 



Table 5. Model Emissivities (vju) band (ergs s 1 sr 1 /H) a for Akari 





model 




N3 




Nl 




S7 




S9W 




Sll 




L15 




L18W 




L24 




N60 


WIDE-S 


WIDE-L 




N160 


u 


flpATT (%1 




3.2/jm 


A.lfim 




7.0/tm 




9.0/im 


ll.OjUm 




15. /im 




18. /im 




24. (im 




65. /jm 




80./im 


140. /im 


160. fim 


0.5 


0.47 


3 


,75e-28 


1 


68e-28 


4. 


,37e-27 


.3 


,55e-27 


3 


49e-27 


5 


13e-27 


5. 


48e-27 


5 


,88e-27 


4 


,46e-26 


8 


.43e-26 


1.32e-25 


1 


26e-25 


0.5 


1.12 


1 


03e-27 


4. 


67e-28 


1 


03e-26 


7 


,54e-27 


6 


,36e-27 


6 


,71c-27 


6. 


38e-27 


6 


,18e-27 


4 


,20e-26 


8 


.02e-26 


1.28e-25 


1 


22e-25 


0.5 


1.77 


1 


69e-27 


7. 


70e-28 


1 


62e-26 


1 


.16e-26 


9 


26e-27 


8. 


,30e-27 


7. 


28e-27 


6 


,47e-27 


.3 


.99e-26 


7 


.73e-26 


1.25e-25 


1 


20e-25 


0.5 


3.19 


3 


,16c-27 


1 


,46e-27 


2 


,95e-26 


2 


.09e-26 


1 


,63e-26 


1 


,29e-26 


1. 


06e-26 


8 


,57e-27 


3 


.62e-26 


7 


,04e-26 


1.18e-25 


1 


,13e-25 


0.5 


4.58 


i 


,68e-27 


2 


,15e-27 


4. 


33e-26 


.3 


.05e-26 


2 


36e-26 


1 


,79e-26 


1 


,42e-26 


1 


.10e-26 


3 


.25e-26 


6 


.31c-26 


1.10e-25 


1 


,06e-25 


1.0 


0.47 


7 


.50e-28 


3 


35e-28 


8. 


,75e-27 


7 


,12e-27 


7 


.Olc-27 


1 


,04e-26 


1 


,12e-26 


1 


21e-26 


1 


.23e-25 


2 


.03e-25 


2.48e-25 


2 


25e-25 


1.0 


1.12 


2 


,05c-27 


9. 


,33e-28 


2 


06e-26 


1 


51e-26 


1 


,27e-26 


1 


35e-26 


1. 


,30e-26 


1 


.27e-26 


1 


.16e-25 


1 


.93e-25 


2.40e-25 


2 


,19e-25 


1.0 


1.77 


3 


,37c-27 


1 


,54e-27 


3. 


25e-26 


2 


.31e-26 


1 


,85e-26 


1 


67e-26 


1. 


,47e-26 


1 


.32e-26 


1 


.llc-25 


1 


.87e-25 


2.36e-25 


2 


,16e-25 


1.0 


3.19 


6 


.32e-27 


2 


91e-27 


5. 


90e-26 


4 


.18e-26 


3 


26e-26 


2 


,59e-26 


2 


13e-26 


1 


.74e-26 


1 


01e-25 


1 


.71e-25 


2.22e-25 


2 


04e-25 


1.0 


4.58 


9 


,34e-27 


4. 


30e-27 


8. 


66e-26 


6 


.10e-26 


1 


72e-26 


3 


60e-26 


2. 


86e-26 


2 


.23e-26 


8 


97e-26 


1 


.55e-25 


2.09e-25 


1 


93e-25 


3.0 


0.47 


2 


,25c-27 


1 


01e-27 


2 


63e-26 


2 


.14e-26 


2 


,12e-26 


3 


,24e-26 


3 


56e-26 


3 


.93e-26 


5 


.66e-25 


7 


.27e-25 


5.93e-25 


5 


.01e-25 


3.0 


1.12 


6 


,16e-27 


2 


80e-27 


6. 


17c-26 


4 


.54e-26 


3 


84e-26 


4 


,18e-26 


4. 


08e-26 


4 


.08e-26 


5 


.38e-25 


6 


.98e-25 


5.78e-25 


1 


89e-25 


3.0 


1.77 


1 


01e-26 


4. 


62e-27 


9. 


76c-26 


6 


.95e-26 


5 


58e-26 


5 


,12e-26 


4. 


60e-26 


4 


.23e-26 


5 


.19e-25 


6 


.80e-25 


5.71e-25 


4 


,84e-25 


3.0 


3.19 


1 


.89e-26 


8. 


,73e-27 


1 


.77e-25 


1 


.25e-25 


9 


,80e-26 


7. 


,88e-26 


6. 


56e-26 


5 


.48e-26 


4 


,73c-25 


6 


,28e-25 


5.43e-25 


4 


.61e-25 


3.0 


4.58 


2 


,80e-26 


1 


29e-26 


2. 


60e-25 


1 


.83e-25 


1 


,42e-25 


1 


,09e-25 


8. 


,77e-26 


6 


.95e-26 


4 


.24e-25 


5 


.76e-25 


5.17e-25 


1 


42e-25 


5.0 


0.47 


.3 


,75e-27 


1 


68e-27 


4. 


39e-26 


.3 


.59e-26 


3 


57c-26 


5. 


,58e-26 


6. 


20e-26 


6 


.93e-26 


1 


.09e-24 


1 


.25e-24 


8.52e-25 


6 


96e-25 


5.0 


1.12 


1 


.03e-26 


4. 


,67e-27 


1. 


03e-25 


7 


.57e-26 


6 


,43e-26 


7 


.12e-26 


7 


05e-26 


7 


.16e-26 


1 


,04e-24 


1 


,20e-24 


8.32e-25 


6 


.80e-25 


5.0 


1.77 


1 


.69e-26 


7. 


,70e-27 


1. 


63e-25 


1 


.16e-25 


9 


33e-26 


8. 


,68e-26 


7 


90e-26 


7 


.37e-26 


1 


OOc-24 


1 


.17e-24 


8.23e-25 


6 


.74e-25 


5.0 


3.19 


3 


,16e-26 


1 


45e-26 


2. 


95e-25 


2 


.09e-25 


1 


,64e-25 


1 


33e-25 


1 


12e-25 


9 


.45e-26 


9 


.18e-25 


1 


.09e-24 


7.85e-25 


6 


,45e-25 


5.0 


4.58 


4 


67e-26 


2 


,15e-26 


4. 


33e-25 


3 


,05e-25 


2 


37e-25 


1 


,83e-25 


1. 


,48e-25 


1 


.19e-25 


8 


.29e-25 


1 


.01e-24 


7.53e-25 


6 


.22e-25 


10.0 


0.47 


7 


50e-27 


3 


35e-27 


8. 


82e-26 


7 


.22e-26 


7. 


,25e-26 


1 


,18e-25 


1. 


34e-25 


1 


.52e-25 


2 


.48e-24 


2 


.50e-24 


1.36e-24 


1 


,06e-24 


10.0 


1.12 


2 


.05e-26 


9 


33e-27 


2. 


06e-25 


1 


.52e-25 


1 


,30e-25 


1 


,48e-25 


1 


,50e-25 


1 


,55e-25 


2 


37e-24 


2 


.41e-24 


1.33e-24 


1 


04e-24 


10.0 


1.77 


3 


.37e-26 


1. 


54e-26 


3. 


,25e-25 


2 


.32e-25 


1 


,87e-25 


1 


,79c-25 


1 


,66e-25 


1 


.59e-25 


2 


.30e-24 


2 


.36e-24 


1.32e-24 


1 


04e-24 


10.0 


3.19 


6 


.32e-26 


2 


91e-26 


5. 


,90e-25 


4 


.18e-25 


3 


,28e-25 


2 


,71c-25 


2. 


31e-25 


2 


.00e-25 


2 


,12c-24 


2 


.21e-24 


1.26e-24 


9 


,96e-25 


10.0 


4.58 


9 


,34e-26 


4. 


30e-26 


8. 


,67e-25 


6 


.llc-25 


4 


,74c-25 


3 


,72c-25 


3 


05e-25 


2 


.50e-25 


1 


93e-24 


2 


.06e-24 


1.22e-24 


9 


65e-25 


100.0 


0.47 


7 


51e-26 


3 


,37e-26 


9. 


26e-25 


7. 


97e-25 


8. 


,76c-25 


1 


,72c-24 


2. 


lie- 24 


2 


63e-24 


2 


55e-23 


1 


.77e-23 


4.84e-24 


3 


,36e-24 


100.0 


1.12 


2 


,05e-25 


9 


37e-26 


2 


.lie- 24 


1 


.59e-24 


1 


,44e-24 


1 


.99e-24 


2. 


,22e-24 


2 


60e-24 


2 


47e-23 


1 


.72e-23 


4.75e-24 


3 


,30e-24 


100.0 


1.77 


.3 


.38e-25 


1 


,54e-25 


3 


30e-24 


2 


,39e-24 


2 


OOc-24 


2 


25e-24 


2 


33e-24 


2 


,57e-24 


2 


.43e-23 


1 


.70e-23 


4.73e-24 


.3 


,29e-24 


100.0 


3.19 


6 


33e-25 


2. 


92e-25 


5. 


96e-24 


4 


,26e-24 


3 


,41e-24 


3 


16e-24 


2 


95e-24 


2 


91e-24 


2 


.29e-23 


1 


.62e-23 


4.60e-24 


.3 


,20e-24 


100.0 


4.58 


9 


36e-25 


1 


31e-25 


8. 


,74e-24 


6 


,19e-24 


4. 


,88e-24 


4 


20e-24 


3. 


69e-24 


.3 


,36e-24 


2 


.15e-23 


1 


,55e-23 


4.51e-24 


3 


,15e-24 



a Additional models, with spectra, are available online at http://www.astro.princeton.edu/~draine/dust/irem.html 
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Table 6. Model Emissivities (vj u )ba,nd (ergs s 1 sr 1 /H) a for Herschel 



model PACS 1 PACS 2 PACS 3 SPIRE 1 SPIRE 2 SPIRE 3 

U ?pah (%) 75pm 110/im 170/xm 250/^m 360/im 520/im 



0.5 


o 


47 


4 


21e- 


26 


7 


92e- 


26 


9 


61e- 


26 


5 


71e- 


26 


2 


05e- 


26 


5 


80e- 


27 


0.5 


1 


12 


3 


97e- 


26 


7 


57e- 


26 


g 


31e- 


26 


5 


58e- 


26 


2 


Ole- 


26 


5 


70e- 


27 


0.5 


1 


77 


3 


78e- 


26 


7 


34e- 


26 


9 


17e- 


26 


5 


53e- 


26 


2 


OOe- 


26 


5 


68e- 


27 


0.5 


3 


19 


3 


43e- 


26 


fj 


71e- 


26 




61e- 


26 


5 


28e- 


26 


1 


92e- 


26 


5 


50e- 


27 


0.5 


4 


58 


3 


05e- 


26 


ft. 


05e- 


26 




09e- 


26 


5 


07e- 


26 


1 


87e- 


26 


5 


39e- 


27 


1.0 


o 


47 


1 


13e- 


25 


1 


78e- 


25 


1 


73e- 


25 




74e- 


26 


2 


83e- 


26 


7 


50e- 


27 


1.0 


1 


12 


1 


08e- 


25 


1 


71e- 


25 


1 


68e- 


25 




55e- 


26 


2 


77e- 


26 


7 


38e- 


27 


1.0 


1 


77 


1 


03e- 


25 


1 


67e- 


25 


1 


I ■- 


25 




48e- 


26 


2 


76e- 


26 


7 


35e- 


27 


1.0 


3 


19 


9 


33e- 


26 


1 


53e- 


25 


1 


57e- 


25 


8 


13e- 


26 


2 


66e- 


26 


7 


13e- 


27 


1 n 


1 


oo 


o 
o 


oze- 


Of; 
ZO 


i 

i 


4ue- 


or. 
Zo 


i 

i 


4toC- 




7 


ooe- 


Oft 

zo 


o 
z 


Due- 


Oft 

zo 


i 


uie- 


07 
Z i 


3.0 





47 


4 


96e- 


25 


5 


69e- 


25 


3 


94e- 


25 


1 
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a Additional models, with spectra, are available online at http://www.astro.princeton.edu/~draine/dust/irem.html 
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Fig. 16. — IRAC color-color plots for dust emission. Symbols with error bars: diffuse emission extracted by Flagey 
et al. (2006) for 5 GFLS fields and for one GLIMPSE field. Squares: present models for Milky Way dust with 
different values of gpAH = 0.47%, 1.1%, 1.8%, and 4.6%, for starlight intensities U < 10 3 so that single-photon 
heating is dominant. The dust model appears consistent with the observed diffuse emission in IRAC 3.6, 5.7, and 
7.9pm bands, but observed emission in the 4.5 /im band appears to be about 50% stronger than given by the model 
emissivities. 



Infrared MidPlane Survey Extraordinaire (GLIMPSE). The model emission spectra, convolved with 
the response function for the IRAC bands, can be compared to the observed colors of the diffuse 
emission. 

Figure PT6a shows the colors of the observed emission in the IRAC 3.6 /xm, 5.7 /xm, and 7.9 /im 
bands. The filled square show the color calculated for the present dust model with the value of 
(/pah ~ 4.6% that appears t o be applicab l e to t he dust in the Milky Way and other spiral galaxies 
with near-solar metallicity (IDraine et al.1 120061 ) . The model with (/pah ~ 4.6% appears to be in 
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Fig. 17. — Upper panel: estimator (l22jl for total infrared flux (TIR) divided by actual TIR. The estimator (f22|) is 
within ±10% of the actual TIR (within the dotted lines) for 0.1 < U < 10 2 and 10 4 < U < 10 6 . Lower panel: the 
TIR estimator from Dale & Helou (2002), applied to the present dust models. 



good agreement with the Milky Way observations in Figure [T6a. First of all, (^i^)s. 7/(^^)7. 9 is 
within the range of observed values, indicating that the adopted C a b s (A) has approximately the 
correct shape over the 5.5-8.5 /im wavelength range. Second, the ratio of 3.6 /jm emission divided 
by 8.0 /im emission is close to the observed value. The emission into the IRAC 3.6 ^m band is 
sensitive to the abundances of the very smallest PAHs, with Ac < 75 (see Figure [6]), we can 
therefore conclude that the size distribution adopted by the present model has approximately the 
correct amount of PAH mass in the interval 25 < Ac < 75, at least relative to the PAH mass in the 
interval 25 < Ac < 10 3 that radiates efficiently into the IRAC 8 /um band. The ratio of the 3.6 fim 
emission to the 8 fim emission of course is sensitive to the assumed PAH charge state (see Fig. [6]) . 
The present model assumes a mixture of ionization conditions that is intended to be representative 
of the local interstellar medium (see Figure [8]) . If the actual ionized fraction for the PAHs is higher 
than in our model, then the actual mass fraction in the 25 < Ac < 75 would have to be even higher 
than the adopted size distribution (shown in Figure [IT]) . Conversely, if the PAH neutral fraction is 
actually higher than in our model, then the PAH mass fraction in the 25 < Ac < 75 range could 
be reduced. 

Observations of external galaxies by ISO and Spitzer Space Telesc ope have been interpreted as 



show ing diffuse nonstellar emission in the 2.5-5 /jm wavelength range (|Lu et al.l 120031 ; iHelou et al 
20041 ) . consistent with thermal emission from dust with Q a b s oc u 2 and temperatures in the ~ 750- 



1000 K range. The GFLS and GLIMPSE images also appear to contain diffuse emission in the 
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4.5 //m band (Flagey et al. 2006). Figure [16b shows the observed ratio of 4.5/xm emission relative 
to 7.9/im emission as well as the model values. The observed 4.5 fim emission is perhaps 50% 
stronger (relative to the 8.0 fim emission) than the models. 

In many cases it is desirable to try to estimate the total infrared luminosity from the observed 
fluxes in the IRAC and MIPS bands. The following weighted sum of IRAC 7.9/im and the 3 MIPS 
bands is a reasonably accurate estimator for -Ltir : 



L TIR « 0.95(^) 7 . 9 + 1.15(1/^)24 + (vL v ) n + (vL„) 



160 



(22) 



Figure [17] plots this estimator relative to the actual Ltir for our dust models, showing that eq. 
(122H allows the total infrared luminosity to be estimated to within ~ 10% for our dust models 
heated by starlight with 0.1 < U < 10 2 or with 10 4 < U < 10 6 . The accuracy becomes somewhat 
lower for 10 2 < U < 10 4 (because of the factor of 3 gap in wavelength between the 71 fim band 
and the 24 /im band), but the worst-case error is only ~ 30%, occuring for U ~ 10 3 where the FIR 
emission peak falls at ~ 40//m, halfway between the MIPS 24/im and 71/xm filters, with the result 
that MIPS photometr y underestimates the actual power. Also plotted in Figure [T7] is the TIR 
esti mator proposed by iDale &: Heloul (|2002l ) using only MIPS photometry, applied to our models. 
The lDale fc Heloul ( 20021 ) luminosity estimate is accurate to within ±25% for U < 10 2 . 



8. Emission Spectra for Dust Models: Distribution of Starlight Intensities 



The dust grains in a galaxy will be exposed to a wide range of starlight intensities. The bulk 
of the dust in the diffuse interstellar medium will be heated by a general diffuse radiation field 
contributed by many stars. However, some dust grains will happen to be located in regions close 
to luminous stars, such as photodissociation regions near OB stars, where the starlight heating the 
dust will be much more intense than the diffuse starlight illuminating the bulk of the grains. 

In principle, one could construct a model for the distribution of stars and dust in the galaxy, 
and solve the radiative transfer problem to determine the heating rate for each dust grain. This, 
however, requires many uncertain assu mptions, as well as heavy numerica l calcu l ations to solve 



the ra d iative trans f er problem [see, e.g. , IWitt &: Gordon! (jl996l ): ISilva et al.1 (|1998l ): iPopescu et al 



d2OO0l ): lTuffs et alJ (|2004h : IPiovan et all feood )]. rlere we take a much simpler approach, and assume 
a simple parametric form for the fraction of the dust mass exposed to a distribution of starlight 
intensities U described by a delta function and a power-law distribution for U m ™ < U < C/ max : 

dM dust n _ A)l/r ^ TT TT , , («-l) 



dU 



(1 - -f)M dust 5(U - U min ) + 7 M dust 



l-a 
max 



-u- 



(23) 



where dM dnst is the mass of dust heated by starlight intensities in [U, U + dU]; M dust is the total 
mass of dust, (1 — 7) is the fraction of the dust mass that is exposed to starlight intensity [/ mm , 
and a is a power -law index. This function al form is similar to the power-law distribution used by 
Dale et al.1 (|200ll ) and lDale &: Heloul ([20020 . except that we have added a delta function component 
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that contains most of the dust mass; the delta- function term is intended to represent the dust in 
the general diffuse interstellar medium. 




Fig. 18. — Emission spectra for models with MW dust with (a) (/pah = 4.6%, (7mm = 1, (b) §pah = 4.6%, 
LUn = 10, (c) (jpah = 1.77%, U min = 1, (d) (jpah = 1.77%, Umin = 10. All models with [/max = 10 6 , a = 2. Values 
of 7 are indicated. 



The star l ight d istribution function eq. (f23|) has parameters U m \ n , t/ max , a, and 7. However, 



Draine et al.l (j2006l ) find that the SEDs of galaxies in the SINGS survey appear to be satisfactorily 
reproduced with fixed a = 2 and U max = 10 6 , and we will therefore use these fixed values for 
a and U max . For a = 2 the fraction 7 of the dust mass that is exposed to starlight intensities 
t^min < U < t/max has equal amounts of infrared power per unit log U. Each of the grain models in 
Table [3] has a different value of qpaHj so the parameter ^pah serves as a proxy for the dust model. 
This leaves us with four free parameters: -Md us t, (?pah> ^mim 7- The shape of the dust emission 
spectrum is determined by only 3 free parameters: qpah, CAnin, and 7. Examples of emission spectra 
are shown in Figure [181 
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Fig. 19. — Ratios of (vL v )7.g/LTm and (v L v )24 / Ltik for dust models with (/pah = 0.47% and 4.6%, plotted against 
the ratio (vL v )n / (vL v )\u) ■ The model grids are labelled by the minimum starlight intensity (7min and the fraction 
7 of the dust exposed to starlight with U > U m in (see text). 



9. Estimating q P AU, U min , 7, /pdr, and M dust 

Photometry obtained with the IRAC and MIPS cameras on Spitzer Space Telescope can be 
used to estimate the free parameters qpAH, U m i n , and 7, the fraction /pdr of the dust infrared 
emission radiated by dust grains in photodissociation regions where U > 10 2 , and the total dust 
mass Md U st- The flux measured in the IRAC 3.6/xm band is almost entirely due to starlight, and 
therefore can be used to remove the starlight contribution to the 7.9/mi and 24/im bands, with the 
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nonstellar flux densities estimated to be 



F-(7.9 M m) 
F^(24/mi) 



F v (7.9fim) - 0.232^(3.6/^) 
F„(24/fln) - 0.032^,(3.6/^1) 



(24) 
(25) 



where the coefficients 0.232 and 0.032 are from Helou et al. (2004). The stellar contribution to the 
71 /jm and 160 /im bands is negligible. 

For a given dust type (i.e., a given value of (/pah) there is a two dimensional family of emission 
models, parameterized by U m { n and 7. Figures [T9a and b show the flux into the MIPS 24/mi 
band normalized by the total dust luminosity Ltir ; plotted against the ratio of the fluxes into 
the MIPS 71//m and 160/mi bands. The two-dimensional family of emission models is shown, for 
0.1 < t/min < 20, and < 7 < 0.20. The models with 7 = (the lower boundary of the model grid) 
are models with dust heated by a single starlight intensity: (U) = U m i n . 

For the 7 = models, the ratio (^-^^/-^tir is essentially independent of U m \ n for U m i n < 2, as 
already seen in Figure [T5l the emission into the 24/zm band is dominated by single-photon heating of 
small grains. For 7 > 0, some fraction of the dust is exposed to starlight intensities U > 10, and the 
fraction of the total power that enters the 24/im band increases. For the particular size distributions 
that have been adopted for these models (see Figure [TTj) . the models with <7pah = 0.47% and 4.6% 
(Figures [T9"a.b differ by only a factor ~ 1.5 in (^L^)24/-^tir in the single-photon heating limit 



The emission into the 7.9/im band is a different story. Figures [T9b.d show that in the single- 
photon heating limit, the fraction of the power radiated into the 7.9/iin band is essentially pro- 
portional to (/pah: as (/pah increases from 0.47% to 4.6%, (^L^^g/Lpm increases from 0.02 to 
0.2. As 7 is increased, some additional 7.9/mi emission is produced; this results in a noticeable 
fractional increase for the (/pah = 0.47% model, but for the (/pah = 4.6% model the additional 
emission produces a very small fractional change in {vL v )-j$/ L^ir, which remains dominated by 
single-photon heating. 

To determine qpKR, U m \ n , 7, and Md ust , the best procedure is to vary all the parameters to 
find the dust model which comes closest to reproducing the photometry. This procedure is used by 
Draine et al. (2006) in their study of the dust properties of the SINGS galaxy sample. 

However, it is also possible to use a graphical procedure to find values of (/pah, ^mim and 7 
that are consistent with the observed data. We assume that we have Spitzer photometry in 5 bands: 
the 3.6/mi and 7.9^m bands of IRAC, and the 24, 71, and 160/mi bands of MIPS. The 3.6/mi band 
observations are used with eq. ([24"]) and ([25]) to obtain (i ? " s )7.g and (i^ s )24- We define three ratios 
of observables: 



(tfmin < 1, 7 = 0). 
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Rn ^ . (28) 

For starlight intensities 0.1 < U < 10 2 , the bulk of the power radiated by dust emerges in the 
50-200 /iin wavelength range, and therefore the total dust luminosity is approximately proportional 
to {{vF u )i\ + (vF v )i§q\. The ratio Rf\ is sensitive to the temperature of the a > 0.01/iin grains 
that dominate the far-infrared emission; P71 is thefore an indicator for the intensity of the starlight 
heating the dust. 



9.1. Determining the PAH Fraction (/pah 

P7.9 is proportional to the fraction of the dust power radiated in the PAH features. Because 
the 7.9/im emission is almost entirely the result of single-photon heating, P7.9 depends very weakly 
on the starlight intensity. 

Figure [20] shows P7.9 vs. P71 for a sequence of grain models, each with a different PAH 
abundance qpah- The model with the lowest value of (/pah naturally has very low values of P7.9. 
Because of the very low PAH abundance in Figure [20b. a very small amount of dust exposed to high 
values of U can enhance the 8 fim emission, and therefore P7.9 increases notably as 7 increases from 
0, with increasing amounts of dust exposed to radiation intensities between U m \ n and U max = 10 6 . 
However, as qpah is increased, P7.9 becomes less sensitive to 7, because single-photon heating of 
the PAHs produces 7.9/im emission even when U is small. 

Based on modeling many spectra, we expect that the 7.9/im emission will primarily arise from 
single-photon heating, i.e., P7.9 will be close to the value calculated for 7 = 0. The approximate 
value of (/pah can therefore be determined by finding a dust model among Figures [20k- g for which 
the observed (P^ijP^.g) point falls just above the 7 = curve. The location (^71^7.9) on the 
appropriate plot also gives an estimate for U m i n and 7, but these parameters are better determined 
using the 24/iin emission, as described below. 



9.2. Determining U m i n and 7 

P24 is proportional to the fraction of the dust power radiated near 24 /iin. As seen in Figure 
[7[ 24/im emission can be produced by single-photon heating of a ~ 15 — 40 A grains. For a given 
dust model there is therefore a limiting value of P24 that applies for 0.1 < U < 10 where the bulk 
of the dust power is captured in the 71 and 160 /iin bands. However, when some fraction of the 
grains are heated by starlight with intensities U > 20 (see Figure [T5|) . larger dust grains can be 
heated sufficiently to add to the 24/im emission. As a result P24 is sensitive to the value of 7. 

To estimate the value of 7, one looks for a dust model in Figure [21] with the value of (/pah 
found from Figure 1201 The observed location of (P71, P24) on Fig. [21] then allows U m i n and 7 to be 
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determined. It may happen that (U m i n , 7) estimated from Figure |2T1 differs from (U m i n , 7) indicated 
by Figure [20l When this occurs, it is an indication that a single dust model does not perfectly 
reproduce the observed 7.9/24/71/160/mi colors. We recommend using the values of (J7 m i n) 7) 
estimated from Figure I2T1 as they are less sensitive to the adopted value of ^pah- 
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Fig. 20.— 



IRAC 7.9 urn power relative to MIPS 71 and 160 /tm power versus MIPS 71/MIPS 160 band ratio. 
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Fig. 21. — MIPS band ratios for models. 

9.3. Determining /pdr 

Young stars tend to be in or near to dust clouds. In star-forming galaxies such as the Milky 
Way, we expect that a significant fraction of the starlight emitted by O and B stars will be absorbed 
by dust that happens to be relatively close to the star, so that the starlight intensity is significantly 
above the average starlight intensity. Photodissociation regions (PDRs) are typical examples of 
such environments. 

For the power-law distribution of starlight intensities given by eq. (|23p . we can calculate a 
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Fig. 22. — /pdr, the fraction of the dust infrared luminosity radiated by dust grains in regions 
where U > 10 2 , plotted as a function of the observable (P24 — 0.14JV.g). For the 5 example galaxies, 
the width of the rectangle corresponds to the ±l<r uncertainty in -P24 - O.I4.F7.9. The fitting function 
1.05(^24 — O.I4P7.9 — 0.035) - 75 approximates the model results. 



quantity /pdr> which we define here as the fraction of the total dust luminosity that is radiated by 
dust grains in regions where U > 10 2 : 

/pDR = 7 ln(tWl0 2 ) (29) 

(1 - 7)(1 - U m i n /U max ) + 7ln(^max/^min) 

If U-m ™ and 7 are known (e.g., either determined by finding the best-fit (ffpAH? U-m in , 7) from the 
3-dimensional model space, or by following the graphical procedures described above), /pdr can 
be calculated from eq. (f29j) . Here we show how /pdr can be estimated directly from IRAC and 
MIPS photometry. 

For the Milky Way dust models considered here, PAHs undergoing single-photon heating con- 
vert a small fraction of the absorbed starlight power into 24//m emission, but when a distribution of 
starlight intensities is present, there is additional 24/um emission from high intensity regions. The 
24fim contribution from PAHs depends, of course, on the PAH abundance. We find that for our 
models /pdr can be closely related to the fraction of the power radiated at 24^m after subtraction 
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of the contribution of PAHs to the 24/xm power. We find that the combination (P24 — O.I4.FV.9) 
gives a quantity that is sensitive to U and relatively insensitive to the value of (?pah- 

Figure [22] shows /pdr versus P24 — 0.14P 7 9 for eight different dust models with different values 
°f Qpah and C/ m i n . We see that the curves fall within a narrow band. For these dust models, /pdr 
can be estimated from the observable P24 — O.I4.P7.9 using Figure! 
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Fig. 23. — Value of "J for dust mass estimation (see text). 
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9.4. Examples 

Table 7: Galaxy Examples. 





IC2574 


Mark 33 


NGC 1266 


NGC 3521 


NGC 6822 


D (Mpc) 


4.02 


21.7 


31.3 


9.0 


0.49 


(^>3.e (Jy) a 


0.156 ±0.016 


0.027 ± 0.003 


0.056 ± 0.006 


2.12 ±0.21 


2.20 ±0.22 


{F v ) 7 . (Jy) a 


0.066 ± 0.007 


0.125 ±0.013 


0.087 ± 0.009 


6.23 ±0.62 


1.41 ±0.14 


(F,)24 (Jy) a 


0.278 ±0.013 


0.836 ± 0.034 


0.861 ± 0.035 


5.47 ±0.22 


2.59 ±0.10 


(F v )n (Jy) a 


5.10 ±0.37 


3.82 ± 0.28 


11.45 ±0.80 


57.2 ±4.0 


59.2 ±4.2 


<^>i6o (Jy) a 


10.8 ± 1.4 


3.63 ± 0.46 


9.6 ± 1.2 


217. ±26 


143 ± 17 


P7.9 


.027 


0.196 


0.0435 


0.529 


0.065 


P24 


.0810 


0.460 


0.168 


0.166 


0.061 


Rn 


1.08 


2.33 


2.55 


0.582 


0.944 


Qpah from Fig. [20] 


0.5% 


3.2% 


0.5% 


4.6% 


1.1% 


t/ min from Fig. Hi] 


2.3 


7 


12.5 


1.1 


2 


7 from Fig. [5T] 


0.012 


0.14 


0.029 


0.007 


0.006 


([/) (from eq. [33J 


2.6 


18. 


16. 


1.2 


2.1 


* (g/(ergss- 1 )) 


0.046 


0.061 


.049 


0.054 


0.0465 


log lo (M dust /M ) 


5.79 


6.41 


7.04 


8.10 


5.15 


a Dale et al. (2006) 



To illustrate this, five galaxies are indicated in Figures [20] and ED NGC 1266, NGC 3521, 
NGC 6822, IC 25 74, and Markarian 33. Table [7] gives the measured global fluxes for these galaxies 
(jDale et al,l 120051 ). the derived ratios -P7.9, P24, and Rn, and the values of (/pah, U m i n , and 7 
obtained from the procedure described in § §9. 1 19.21 Below we comment on the individual cases. 

IC 2574: In Figure [20], the very low value of -P7.9 = 0.027 only falls on the model grid in Figure 
[20k . therefore we select (/pah = 0.5%. From Figure [2Th. we find U m - m 2.3 and 7 0.012. 

Mark 33: From Figure 120b . we estimate (/pah ~ 3.2%. The very high value of P24 = 0.46 
requires U m i n ~ 7 and 7 ~ 0.14 on Figure l2l"b in order to have enough hot dust to reproduce the 
observed 24//m flux. The very large value found for 7 suggests that a large fraction of the dust 
heating in Mark 33 is taking place in photodis sociation regions near OB associations, consistent 
with the dwarf starburst nature of this galaxy ([Hunt et al.ll2005l ). 



NGC 1266: The low value of P7.9 ~ 0.044 can only be reproduced in Figure [20] by models with 
small values of gpAH- Of the values of (/pah for which spectra have been computed, (/pah = 0.47% 
is preferred. Turning now to Figure l2~Ta. we find that U m i n ~ 12.5 and 7 ~ 0.029. The values of 
<?PAH; kmin> and 7 found by this graphical e xercise are in good agreement with the best-fit values 
of q PAH = 0.47%, U min = 12, and 7 = 0.029 dDraine et alJbood l. 



NGC 3521: From Figure [20b. we find (/pah ~ 4.6%. Fro m Figure [21k. we obtain U m i n « 1.1 and 
7 ~ 0.007. SCUBA fluxes are also av ailable for this galax y ([Dale et al.ll2005l ): when the 850/mi flux 



is used to constrain the grain model, iDraine et al.l ([2006) find that it is difficult to simultaneously 
reproduce the MIPS and SCUBA photometry. The overall best fit is found for U m i n = 2, 7 = 0.01. 
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NGC 6822: From Figure [20b we estimate qpah = 1.1%. From Figure l2"Tb. we estimate f7 m i n ~ 
2.0 and 7 « 0.006. These are i n agre ement with the best-fit values <jpah = 1-2%, ?7 m i n = 2.0, 
7 = 0.005 found bv braine etaD (|20od ). 



9.5. Estimating the Dust Mass Md ust 

Consider a galaxy at distance D. If is the dust emissivity per H nucleon, then the flux is 

ft -^4 . (30) 

and the dust mass 

M dust =(^)m H ^D 2 . (31) 

If we know the values of {M dvLS t/My) and j v for a given grain model, we can estimate Md us t directly 
from the distance D and the observed flux F v . However, our estimate for j u may depend sensitively 
on our estimates for (/pah> U-mm, and 7. 

A relatively robust approach to dust mass estimation is to note that for a given dust mixture 
the total dust luminosity is proportional to M^ US ^,(U). Define the quantity 

*(«pah,7,CW) = —TT— m KT — . — - , (32) 

where mn is the mass of an H atom, and (?7) is the mean starlight intensity. For the distribution 
function (|23p . with a = 2, (£7) is 

/rr\ _ / UdM dust _ lUmin ln(U max / U min ) 

J "-'"dust 1 ~~ <^min/ L^max 

The dust mass is related to the observed fluxes by 

M dnst = — [(vF v ) 2i + {uF v ) 7l + {uF u ) im ]D 2 . (34) 

For a given dust model, we can calculate ^f, which is shown in Figures 1237 a-g). We observe 
that the values of ^ in Figures l23Ta-g) range only from 0.044 to ~0. 066 g/( ergs s" 1 ), with ^ ~ 
0.055 g/( ergs s _1 ) as a representative value. However, if we already have estimates for <7paH; ^mim 
and 7, it is straightforward to obtain an accurate value for ^(^paH) ^minj7) from Figure | 

To summarize, the procedure for estimating the dust mass M dust is as follows: 



1. Use (F™) 7 ,g, (i ? " s )24, (F v )n, and {F u ) wo to calculate P7.9, P 2 4, and R 



71- 



2. Use R71 and P7.9 to estimate qpah by finding the value of gpAH such that (-R71, P7.9) falls just 
above the 7 = curve in Figure [20j 
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3. Using this value of (/pah? locate (^71,^24) on Figure [2T1 to determine U m \ n and 7. 

4. Using the values of (/pah, f^min, and 7, use Figure [23] to find 

5. Use U m \ n and 7 (and C/ max = 10 6 ) to calculate (U) from eq. ([33]) . 

6. Calculate M^ ust using eq. ([MP - 

We have estimated (£/}, Vl/, and the resulting dust mass -Md ust for the 5 sample galaxies, with 
results given in Tabled 



10. Discussion 



10.1. Origin of near-IR Continuum Opacity of PAHs 



The present models assume that the small-particle end of the grain size distribution is dom- 
inated by PAH particles, as assumed by WD01 and Li & Draine (2001b), and the PAH particles 
are therefore assumed to have a "continuum" component to C a b s in order to be able t o reproduce 



2-6 /i m emission seen by ISO (|Lu et al.ll2003l ) and emission in the IRAC 4.5 /im filter (jHelou et al 



2004). Continuum emission also appears to be present at other wavelengths, e.g., between the 12.6 
and 16.45 /im features. In the 2-6 p wavelength range, even very small carbonaceous grains are 
arbitrarily taken to have C a b s equal to 1% of the continuum that would be produced by free electron 
absorption in graphite, added to the discrete absorption bands expected for PAHs. The nature of 
this conti nuum absorption is unclear. Perhaps some PAHs, e.g., tubular PAHs with appropriate 
chirality (|Zhou et al.ll2006l ) have zero bandgap as in graphite. 



We are assuming that the PAHs are responsible for the 2-5/iin continuum. lAn Sellgren 



(|2003l ) find that the 2/im continuum emission and 3.29/iin PAH emission have different spatial 
distributions in NGC 7023. Perhaps the 2/im emission is strongest in regions where the PAHs have 
been dehydrogenated, thereby suppressing the 3.29/im emission. 

The assumption that the continuum emission originates in PAHs is consistent with the fact 
that the 9.8 /im and 18 /jm silicate features are generally not seen in emission, except from regions 
(dusty winds, compact H II regions, dust near the Trapezium stars in Orion) where the starlight 
i s thoug ht to be sufficiently intense to heat a ~ 0.1/iin grains to T > 200 K. However, iLi Draine 
(|2001bl ) showed that the observed emission from the diffuse interstellar medium in the Milky Way 
would in fact allow larger amounts of ultrasmall silicate grains than had been estimated previously, 
because the 9.7 /im silicate emission feature can be "hidden" between the 8.6 /im and 11.3 /im PAH 
features. 



In summary, the origin of the 2-5/im continuum emission remains uncertain. 
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10.2. Near-IR Absorption of PAH Ions: Astrophysical Implications 



The far-red to near-IR opacity of PAH ions described in iLi fc Draind (2001a) was d erived 
from earlier experimental data for a small number of small PAH ions [see ISalama et al.l (11 996) 



Mattioda et al 



and ref erences therein]. More recent laboratory measurements were carried out by 
(|2005bl ) for 27 PAH cations and anions ranging in size from C14H10 to CsoH22- The newly-measured 
absorption of PAH ions at 0.77-2 .5 urn is considerabl y higher than previously reported. This further 
supports our earlier conclusion (|Li &i Draind l2002bl ) that visible and near-IR photons are able to 
excite PAHs to temperatures high enough to emit the mid-IR band s, in agreement with the detection 



of th e PAH features in interstellar regions lackin g UV photons (lUchida et al.1 Il998l : iPagani et al 



1993)0 and in agreement with the observation bv lSellgren et al.l ()l990h that the ratio of the IRAS 
12 fim emission (to which the PAH features are the dominant contributor) to the total far-IR surface 
brightness is independent of the effective temperature T e g of the exciting stars, for 24 reflection 
nebulae for 5000 K < T cS < 33000 K. 

This also relates to our understanding of the origin of interstellar PAHs. At present, the origin 
and evolution of interstellar PAHs are not very clear. Suggested sources for interstellar P AHs include 
the injection (into the interstellar medium) of PAHs formed in carbon star outflows (|Juralll987l : 
Latter! 1 199 ll ). However, the PAH emission features are commonly not seen in the mid-IR spectra 
of C-rich AGB stars. Th e few C stars which di s play the PAH featur es all have a hot companion 
which emits UV photons ([Speck &: Barlowlll997l ; iBoersma et al.l l2006l). It has been suggested that 
PAHs are present in all C stars but they are simply not excited sufficientl y to emit at mid-IR due 
to lack of UV photons. The new visible-IR cross sections for PAH ions iMattioda et al.l ((2005b) 
suggest that PAHs can be excited even by the starlight from a C star, in which case the absence of 
PAH emission places a limit on the abundance of small PAHs in these outflowsl^l 



10.3. Detectability of 6.2/mi Absorption by Interstellar PAHs 

The dust model appropriate to Milky Way dust is thought to have qpar ~ 4.6% of the dust 
mass in PAHs with Ac < 10 3 C atoms. The carbonaceous grains with a < 50 A contain 52 ppm 



11 PAH emission features have a lso been detected in UV- poor dust debris disks around main-se quence stars 
SAO 20 6462 [Sp=F8V, T cS » 6250 K; ICoulson fc Waltherl <| 19951 )] and HD 34700 [Sp=G0V, T cH « 5940K: ISmith et al ' 



(|2004bl )l However, i n an e xtensive Spitzer IRS spectroscopic survey of 111 T Tauri stars in the Taurus star- forming 



region, 



Furlan et al. ( 20061 ) found that the PAH emission bands are not seen in dust disks around T Tauri stars of 



spectral type later than Gl. 

Jura et al. (2006) reported detection of PAH emissio n features in HD 233517, an evolved oxygen-rich K2III red 



giant (T e g w 4390 K) with circumstellar dust. lJura (2003) argue d th at the IR excess around HD 233517 is unlikely to 
be produced by a recent outflow in a stellar wind. lJura et al. [ |2006t ) hypothesized that there is a passive, flared disk 



orbiting HD 233517 and the PAH molecules in the orbiting disk may be synthesized in situ as well as having been 
incorporated from the ISM. 
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of C per H nucleoli, divided approximately equally between neutral and ionized PAHs. For this 
model, we predict a 6.22^/m absorption feature with integrated strength 

/ At dX- 1 = 6.7 x 10- 23 ( 9m£L) cm/H . (35) 



Chiar & Tielend (j200ll ) found an upper limit / At dX' 1 < 0.8 cm" 1 toward Cyg OB2 No. 12. For 
the estimated Ah = 1.9 x 10 22 cm -2 , this gives an upper limit 

J At dX~ l < 4.2 x 1(T 23 cm/H , (36) 

somewhat smaller than our estimate (1351). 



While the present dust model has 52 ppm C in PAHs with Nq < 5 x 10 4 , the 6.22//m emission 
is produced mainly by PAHs with Ac < 500, which account for only 35 pp m. The 6.2/xm absorption 



contributed by PAHs with Ac < 500 in our model is consistent with the IChiar &: Tielend (|200ll ) 
upper limit. 

The absolute PAH band strengths in Table [U are only estimates. Because the 3-25/iin emission 
spectrum for U = 1 is almost entirely the result of single-photon heating, the emission spectrum 
would be almost unaffected if all the PAH band strengths were uniformly reduced by a common 
factor. However, based on the comparison with theoretical calculations shown in FigureEl the actual 
band strengths should be similar in magnitude to the values adopted in Table [Q Spectroscopy of 
heavily extinguished stars using the IRS instrument on Spitzer Space Telescope should be carried 
out to attempt to detect the predicted 6.22/iin absorption from interstellar PAHs. 



10.4. Deuterated PAHs 



Since the aromatic C-D bond zero-point energy is about 30% lower than that of the C-H bond, 
H loss is favored over D loss. The PAH D/H ratio might therefore be substantially larger than the 
total D/H ~ 2 x 10 -5 ( Allamandola et al.lll989bl ). The observed sightline-to-sightline variations 
of (D/H) in interst ell ar gas h a ve be e n inte rpreted as requiring depletion of D by interstellar dust 



(jLinskv et al.l 120061 ) . iDraind (|2004l . 120061 ) has suggested that this depletion of D may take via 



i ntera ctions of PAH ions with the gas, and may imply D/H~ 0.3 in interstellar PAHs. iPeeters et al 



([2004]) report tentative detection of emission at 4.4 and 4.65^m from deuterated P AHs from the 



Orion Bar and M17 PDRs. If the emission is in fact due to C-D stretching modes, IPeeters et al 
(|2004h estimate D/H « 0.17 ±0.03 for the emitting PAHs in the Orion Bar and 0.36 ±0.08 in M17. 



Our model for the PAH absorption cross section is based on observed emission spectra, and 
some of the modeled absorption may be associated with C-D bending modes if interstellar PAHs 
are appreciably deuterated. Because the observations of 4.4 and 4.65//m emission in the Orion Bar 
and M17 are tentative, we do not attempt here to include any absorption and emission that might 
result from C-D stretching modes in the 4.5/xm region. 
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10.5. Ubiquity and Absence of PAHs in Astrophysical Regions: Rationale for 

Variable ^pah 

The ISO and Spitzer imaging and spectroscopy have revealed that PAHs are a ubiquitous 
feature of both the Milky Way and external galaxies. Recent discoveries include the detection of 
PAH emission in a wide ran ge of systems: dist ant Luminous Infrared Galaxies (LIRGs) with redshift 
z ranging from 0.1 t o 1.2 (|Elbaz et alJ 1200 5); distant Ultraluminous Infrared Galaxies (ULIRGs) 
with redshift z ~2 (|Yan et al.ll2005l ); distant luminous submillimeter galaxies at redshift z ~2.8 



(|Lutz et al.ll2005l ): elliptical galaxies with a hostile environment (containing hot gas of temperature 



10 7 K) where PAHs can be easily destroyed through sp uttering by plasma ions (jKaneda et al 



2005 ); faint tidal dwar f galaxies with metallicit y 
(jlrwin k Madden! bood : lEngelbracht et alJbood l. 



Zq/3 (jHigdon et al.ll2006l ); and galaxy halos 



However, the PAH features are weak or even absent in low-metallicity galaxies and AGN: 



Based on an analy s is of t he mid-IR spectra of 60 galaxies obtained by ground-based observa- 



tions, iRoche et al.l (|199ll ) reported the lack of PAH emission features in AGN. 



The Small Magellanic Cloud (SMC), an irregular dwarf galaxy with a metallicity just ~ 1/10 
of solar, exhibits a local minimum at ~ 12 /im in its IR emission spectrum, suggesting the 



lack of PAHs (|Li Draindl2002d l 



Thuan et al.1 (119991 ) . IPlante & Sauvagd (120021 ). and iHouck et al.1 (120041 ) have shown that the 
PAH features are absent in the ISO and Spitzer spectra of SBS 0335-052, a metal-poor 
(Z ~Z /41) blue compact dwarf galaxy. 



More recently, iHunt et al.1 (|2005l ). Madden et al.1 (J200J), and IWu et al.1 (|2006h performed a 
more systematic investigation of the mid-IR spectra of a large number of low-metallicity 
galaxies (with Z/Zq ranging from 0.02 to 0.6) obtained with ISO and Spitzer. They found 
the PAH features are substantially suppressed in metal-poor dwarf galaxies. 

Using Spitzer IRAC [3.6]-[7.9] colors and Sloa n Digital S k y Sur vey (SDSS) [g-r] colors of 
313 visible-selected SDSS main sample galaxies. iHogg et al.l (|2005l ) found that low-luminosity 
galaxies show a deficiency in PAH emiss ion, with weak evidence for a dependence of the PAH- 
to-stellar radiation ratio on metallicity. iRosenberg et al.l (|2006l ) performed a similar study for 
a statistically complete sample of 19 star-forming dwarf galaxies, but they found a significant 
number of low-luminosity galaxies with very red [3.6]— [7.9] colors, indicating the presence of 
PAHs (and/or hot dust). They also found that the 7.9 fim emission is more strongly correlated 
with the star formation rate than it is with the metallicity. 



Engelbracht et al.l (|2005l ) examined the Spitzer mid-IR colors of 34 galaxies ranging over 2 
orders of magnitude in metallicity. They found that the 8 /xm-to-24 /jm color changes abruptly 
from F„(8/xm)/i^(24/mi)~ 0.7 for galaxies with Z/Z Q > 1/3 to F 1 ,(8 J um)/.F i ,(24/im)~ 0.08 for 
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galaxies with Z/Zq < 1/5-1/3. They attributed this color shift to a decrease in the 7.7/im 
PAH feature at low metallicity. 

The silicate-graphite-PAH model described in this paper has been employed to quantitatively 
determine the fraction qpah of th e dust mass contributed by PAHs for 61 galaxies in the 
SINGS survey (jDraine et al.l 120061 ) . There appears to be a threshold metallicity: galaxies 
with [0/H]= 12 + log 10 (O/H) < 8.1 (i.e., Z < Z /3) have (/pah < 1.5%. 



Knowledge of how ^pah varies as a function of galaxy parameters (such as metallicity, lumi- 
nosity, and morphological type) will have a broad impact - 



It will be essential for properly interpreting the IR c osmological su r veys which now often 
assume a constant PAH strength for galaxies [see, e.g.. lLagache et al.l (200J )]. 



It may provide physical insight into the validity of using the line-to-continuum ratio of the 
7.7 /urn PAH feature as a discriminator between starburst and AGN activity in ULIRGs, i.e 



wheth er the dominant luminosity source of ULIRGs is an AGN or a starburst; see lGenzel et al 
(|l99ah . 

As local analogs of early galaxies at high redshift which must have formed at very low metallic- 
ity, nearby low metallicity systems can provide insight into the early stages of galaxy evolution. 
The low PAH abundances observed in these galaxies reflect the balance of PAH formation 
and destruction, and thereby is a valuable probe of the physical conditions following the onset 
of star formation in primordial gas. 

By comparing with the star-formation rates derived from other tracers, it will allow us to 
quantitatively evaluate the validity of using the IRAC 8 //m photometry as a reliable tracer 
for star-formation rates. It is generally believed that the IRAC 8 /im flux is stronger in 
regions with stronger star-forming activities; however, observations suggest that PAHs are 



destroyed in star-forming regions with very strong an d hard radiation fields (jContursi et al 
2000; iForster Schreiber et alJbood : beirao et allbood ). 



10.6. Diversity of PAH Emission Spectra 

The PAH absorption and emission properties presented in this series of papers are for "astro- 
nomical" PAHs - we fix the peak wavelength and bandwidth of each PAH band to what is observed 
in most astronomical regions. The strength of each PAH band (on a per carbon or hydrogen atom 
basis) is also fixed for a specific charge state to be consistent with what is observed in space, 
measured in the laboratory, or calculated theoretically. 

However, there exist variations in the central wavelength, bandwidth and strength of PAH 
emission bands among different astronomical objects or among different locations within one object: 
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Cohen et al.l (|1989l ) reported that the central wavelength of the 7.7 /im band varies among 
different types of objects: it shifts from shorter wavelengths for objects with heavily processed 
PAHs (e.g. HII regions, reflection nebulae) to longer wavelengths for objects with freshly 
created PAHs (e.g. planetary nebulae). This wavelength shift has been confi rmed by high 



resol ution spectroscopy obtained with the ISO Short Wavelength Spec trometer (IPeeters et al 
2002). Similar shifts have also been seen in other bands [e.g. see Ivan Diedenhoven et al 



(J2004J)] . More recently, iBregman &; Temil (|2005l ) observed a progressive blueshift in the central 
wavelength of the 7.7//m band within 3 reflection nebulae from the edge of the nebulae to 
places closer to the exciting stars. 



Uchida et al.l (|2000l ) reported a progressive broadening of the FWHM of the 7.7//m band 



moving toward the exciting star in the reflecti on nebula vdB 17. T his broadening has also 
been reported for the reflection nebula Ced 21 (jCesarsky et al.ll2000l ). But such a broadening 
has not been seen in other PAH bands. 



• iRoelfsema et al.l (| 19961 ) reported that the 8.6 ^m band of the compact HII region IRAS 18434- 
0242 is almost twice as strong as the 7.7 /im band, while for most objects the 7.7 ^m band is 
stronger than the 8.6 fim band by a factor of ~ 5-10. 



Reach et al.l (|2000i ) reported that the 11.3//m/7.7//m band ratio of the PAH emission spec- 



trum of the quiescent molecular cloud SMC Bl No. 1 in the SMC is much higher than t hat 
of any other Galactic or extragalactic objects, (see Figs. 16 and 17 of lDraine Lil (|200ll )). 



While the band-ratio variations can presumably be accounted for by the combined effects of 
variations in the radiation exciting the PAHs and variations in PAH sizes and ionization state 
(small PAHs emit more strongly at the 6.2, 7.7, 8.6 fim bands, while relatively large PAHs emit 
more strongly at the 11.3 /xm band; neutral PAHs are much stronger 3.3 ^m and 11.3 /im emitters 
and weaker 6.2, 7.7 and 8.6 //m emitters compared to their charged counterparts), we do not intend 
to confront our PAH model with extreme spectral variations, such as the observed spectrum of IRAS 
18434-024 (with the 8. 6 um band almost twice as strong as the 7.7 /im band) which, according to 
Roelfsema et al. (j 19961 ). is possibly rich in non-compact PAHs. Because of the way we design the 
"astronomical" PAH model, the present model is not able to account for the central wavelength 
shifts an d the 7.7 //m band broadening, although experimental investigations have been carried out 



[e.g. see 



Hudgins & Allamandolal (119991 )] 



Hudgins Allamandolal (|2005l ) discuss recent work toward PAH identification, with attention 



to the effects of N substitution. Future, more realistic modeling will requi re more extensive knowl- 
edge o f PAH properties than is now available from laboratory studies [e.g., Hudgins & Allamandola 



2004h : lMulas et al.l iboOfil ) : iMalloci et all hood )}. 



1999); M attioda et al.l (l2005bllar)1 and theoretica l calculations [e.g-. lBakes et al. 



(2001a 



3; 



Malloci et al 
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10.7. Other Grain Models 

The dust model used in this study is, of course, provisional. The model assumes the amorphous 
silicate and carbonaceous material to be in physically separate grains. While this appears to be 
consiste nt with observatio ns, including the observed absence of polarization in the 3.4/xm absorption 



feature (jChiar et al.ll2006l ). it is by no means observationally established. 



Our models have also used the dielectric tensor of crystalline graphite to estimate the infrared 
absorption by carbonaceous grains with a > 0.01/im. Some form of amorphous carbon may provide 
a better approximation for interstellar carbonaceous material that may be formed from agglom- 
eration of PAHs followed by cross-linking due to ultraviolet photolysis and cosmic ray exposure. 
The main effect on the present models would be to do away with the broad absorption feature 
near 30^m that appears in the calculated C a b s (A) for graphitic material, as seen in Figured and 
the corresponding broad emission feature near 30/um when the starlight intensity is high enough to 
heat a > 0.01/um carbonaceous grains to T > 10 2 K. 

In this work we focus on the silicate-graphite-PAH model. For other dust models, if they 
are tailored to reproduce the Milky Way interstellar extinctionand albedo in the optical/UV and 
meanwhile satisfy plausible interstellar abundance constraints, the derived parameters Md us t(U} will 
be nearly the same as derived from the silicate-graphite-PAH model, simply from energy balance 
considerations: the energy emitted in the IR must be absorbed by the same dust in the optical/UV. 
If the dust model reproduces the observed far-infrared emission spectrum of Milky Way dust for 
U ~ 1, then the temperatures of the larger dust grains will scale with changes in U in a manner 
similar to the present model, and therefore the inferred properties of the radiation field (U m i n , 7) 
will presumably be similar. 

Howeve r, the composite model will probably have a lower optical/UV albedo than observed 



(jDwekl 119971 ). Because the extinction per H atom is fixed by observations, the lower albedo will 
result in more total far-infrared emission per H nucleon. Given the uncertainties regarding both the 
intensity of interstellar starlight and the far infrared emission from dust, composite grain models 
must be considered to be viable. Additional observational tests to distinguish between dust models 
are needed. 



11. Summary 

The principal results of this study are as follows: 

1. The wavelength dependence of absorption by interstellar PAHs has been revised to better 
reproduce oberved emission spectra. The adopted C a b s (A) are generally consistent with mea- 
sured or calculated absorption cross sections for PAHs (Figure [2]) . 

2. Temperature distribution functions have been calculated for PAHs of different sizes, heated 
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by starlight with the spectrum estimated for the local interstellar radiation field. These 
temperature distribution functions have been used to compute the time-averaged emission 
spectra for PAHs of different sizes (Figure [5]). 

3. Figures O and [7] show the efficiency with which PAHs radiate absorbed energy into different 
emission features. The strong 7.7 ^m feature is radiated efficiently only by PAHs with Nq < 
10 3 carbon atoms - larger PAHs radiate primarily at longer wavelengths. 

4. Emission spectra have been calculated for dust models that reproduce the average Milky Way 
extinction, using size distributions with different abundances of small PAHs (see Fig. I12j) . 

5. Emission spectra have been calculated for dust mixtures heated by different starlight inten- 
sities (see Fig. [T3"|h 

6. The dust model used in this paper is found to be in good agreement with observed IRAC 3.6/iin/ 
IRAC 7.9/um and IRAC 5.7/iin/ IRAC 7.9/mi band ratios (see Fig. [16]), but underpredicts the 
observed IRAC 4.5/iin/ IRAC 7.9/iin ratio by about a factor 1.5. 

7. The far-infrared and submm emission spectrum calculated for the present model appears to 
be in good agreement with the observed emission from dust in the local Milky Way (see 
Figure [H]). 

8. A prescription (eq. [T7|) is found that allows the total IR emission to be estimated from 
observed fluxes in the IRAC 7.9/mi band and the 3 MIPS bands. For dust heated by starlight 
intensities 0.1 < U < 10 2 , eq. (fT7|) is accurate to within ~ 10%. 

9. Emission spectra have been calculated for dust mixtures heated by power-law distributions 
of starlight intensities (see Fig. I18p . 

10. A distribution function (eq. [23]) is proposed for representing the distribution of starlight 
intensities heating dust in a galaxy, or region within a galaxy. 

11. Graphical procedures are presented that allow dust model parameters to be found using 3 
ratios (-P7.9, and -R71) constructed from fluxes measured in IRAC band 4 and MIPS bands 
1-3. The graphs provided here (Figs. 1201 and I2T1) allow model parameters qpah> kmin> and 7 
to be estimated. 

12. The fraction /pdr of the dust infrared emission that is contributed by dust in photodissocia- 
tion regions with U > 10 2 can be estimated from IRAC and MIPS photometry using Figure 

m 

13. Using a coefficient ^ obtained from Figure [231 the total dust mass Md us t can be estimated 
from fluxes in the 3 MIPS bands using eq. (|34p . 



50 



We are grateful to G. Malloci for making available to us theoretical absorption cross sections 
for a number of PAHs in advance of publication; to E. Peeters for providing SWS spectra of 
the Orion Bar, M17, and NGC 7027; to K. Sellgren for providing IRS spectra of NGC 7023; 
to T. Onaka for providing response functions for Akari; and to A. Poglitsch and M. Griffin for 
providing response functions for PACS and SPIRE. We thank the referee, G. Clayton, for helpful 
comments. BTD was supported in part by NSF grant AST-0406883. AL is supported in part 
by the University of Missouri Summer Research Fellowship, the University of Missouri Research 
Board, and NASA/Spitzer Theory Programs. We are grateful to R.H. Lupton for availability of 
the SM graphics program, used extensively in preparation of this paper. 

REFERENCES 

Adamson, A. J., Whittet, D. C. B., Chrysostomou, A., Hough, J. H., Aitken, D. K., Wright, G. S., 
& Roche, P. F. 1999, ApJ, 512, 224 

Allamandola, L. J., Bregman, J. D., Sandford, S. A., Tielens, A. G. G. M., Witteborn, F. C, 
Wooden, D. H., & Rank, D. 1989a, ApJ, 345, L59 

Allamandola, L. J., Tielens, A. G. G. M., & Barker, J. R. 1985, ApJ, 290, L25 

Allamandola, L. J., Tielens, G. G. M., & Barker, J. R. 1989b, ApJS, 71, 733 

An, J. H. & Sellgren, K. 2003, ApJ, 599, 312 

Asplund, M., Grevesse, N., Sauval, A. J., Allende Prieto, C, & Kiselman, D. 2004, A&A, 417, 751 

Bakes, E. L. O., Tielens, A. G. G. M., & Bauschlicher, Jr., C. W. 2001a, ApJ, 556, 501 

Bakes, E. L. O., Tielens, A. G. G. M., Bauschlicher, Jr., C. W., Hudgins, D. M., & Allamandola, 
L. J. 2001b, ApJ, 560, 261 

Beintema, D. A., van den Ancker, M. E., Molster, F. J., Waters, L. B. F. M., Tielens, A. G. G. M., 
Waelkens, C, de Jong, T., de Graauw, T., Justtanont, K., Yamamura, I., Heras, A., Lahuis, 
F., & Salama, A. 1996, A&A, 315, L369 

Beirao, P., Brandl, B. R., Devost, D., Smith, J. D., Hao, L., & Houck, J. R. 2006, ApJ, 643, LI 

Boersma, C, Hony, S., & Tielens, A. G. G. M. 2006, A&A, 447, 213 

Boulanger, F., Abergel, A., Cesarsky, D., Bernard, J. P., Miville Deschenes, M. A., Verstraete, 
L., & Reach, W. T. 2000, in ESA SP-455: ISO Beyond Point Sources: Studies of Extended 
Infrared Emission, ed. R. J. Laureijs, K. Leech, & M. F. Kessler, 91-98 

Boulanger, F., Boissel, P., Cesarsky, D., & Ryter, C. 1998, A&A, 339, 194 



51 



Bregman, J. & Temi, P. 2005, ApJ, 621, 831 

Cartledge, S. I. B., Lauroesch, J. T., Meyer, D. M., & Sofia, U. J. 2006, ApJ, 641, 327 
Cesarsky, D., Lequeux, J., Ryter, C, & Gerin, M. 2000, A&A, 354, L87 

Chiar, J. E., Adamson, A. J., Whittet, D. C. B., Chrysostomou, A., Hough, J. H., Kerr, T. H., 
Mason, R. E., Poche, P. F., & Wright, G. 2006, ApJ, 651, 268 

Chiar, J. E. & Tielens, A. G. G. M. 2001, ApJ, 550, L207 

Cohen, M., Tielens, A. G. G. M., Bregman, J., Witteborn, F. C, Rank, D. M., Allamandola, L. J., 
Wooden, D., & Jourdain de Muizon, M. 1989, ApJ, 341, 246 

Contursi, A., Lequeux, J., Cesarsky, D., Boulanger, F., Rubio, M., Hanus, M., Sauvage, M., Tran, 
D., Bosnia, A., Madden, S., & Vigroux, L. 2000, A&A, 362, 310 

Coulson, I. M. & Walther, D. M. 1995, MNRAS, 274, 977 

Cunha, K., Hubeny, I., & Lanz, T. 2006, ApJ, 647, L143 

Dale, D. A., Bendo, G. J., Engelbracht, C. W., Gordon, K. D., Regan, M. W., Armus, L., Cannon, 
J. M., Calzetti, D., Draine, B. T., Helou, G., Joseph, R. D., Kennicutt, R. C, Li, A., 
Murphy, E. J., Roussel, H., Walter, F., Hanson, H. M., Hollenbach, D. J., Jarrett, T. H., 
Kewley, L. J., Lamanna, C. A., Leitherer, C, Meyer, M. J., Rieke, G. H., Rieke, M. J., 
Sheth, K., Smith, J. D. T., & Thornley, M. D. 2005, ApJ, 633, 857 

Dale, D. A., Gil de Paz, A., Gordon, K. D., Armus, L., Bendo, G. J., Bianchi, L., Boissier, S., 
Calzetti, D., Engelbracht, C. W., Hanson, H. M., Helou, G., Kennicutt, R., Madore, B. F., 
Martin, D. C, Meyer, M. J., Murphy, E. J., Regan, M. W., Smith, J. D. T., Sosey, M. L., 
Thilker, D. A., & etal. 2006, in preparation, 000, 000 

Dale, D. A. & Helou, G. 2002, ApJ, 576, 159 

Dale, D. A., Helou, G., Contursi, A., Silbermann, N. A., & Kolhatkar, S. 2001, ApJ, 549, 215 
Desert, F.-X., Boulanger, F., k Puget, J. L. 1990, A&A, 237, 215 

Draine, B. T. 1994, in ASP Conf. Ser. 58: The First Symposium on the Infrared Cirrus and Diffuse 
Interstellar Clouds, ed. R. M. Cutri & W. B. Latter, 227-241 

Draine, B. T. 2003, ARA&A, 41, 241 

Draine, B. T. 2004, in Origin and Evolution of the Elements, ed. A. McWilliam & M. Rauch, 
317-335 

Draine, B. T. 2006, in Astronomical Society of the Pacific Conference Series, ed. G. Sonneborn, 
H. W. Moos, & B.-G. Andersson, 58-69 



52 



Draine, B. T., Dale, D., Bendo, G., k et. al. 2006, in preparation, 000, 000 

Draine, B. T. k Lee, H. M. 1984, ApJ, 285, 89 

Draine, B. T. k Li, A. 2001, ApJ, 551, 807 

Duley, W. W. k Williams, D. A. 1981, MNRAS, 196, 269 

Dwek, E. 1997, ApJ, 484, 779 

Elbaz, D., Le Floc'h, E., Dole, H., k Marcillac, D. 2005, A&A, 434, LI 

Engelbracht, C. W., Gordon, K. D., Rieke, G. H., Werner, M. W., Dale, D. A., k Latter, W. B. 
2005, ApJ, 628, L29 

Engelbracht, C. W., Kundurthy, P., Gordon, K. D., Rieke, G. H., Kennicutt, R. C., Smith, J.-D. T., 
Regan, M. W., Makovoz, D., Sosey, M., Draine, B. T., Helou, G., Armus, L., Calzetti, D., 
Meyer, M., Bendo, G. J., Walter, F., Hollenbach, D., Cannon, J. M., Murphy, E. J., Dale, 
D. A., Buckalew, B. A., k Sheth, K. 2006, ApJ, 642, L127 

Finkbeiner, D. P., Davis, M., k Schlegel, D. J. 1999, ApJ, 524, 867 

Flagey, N., Boulanger, F., Verstraete, L., Miville Deschenes, M. A., Noriega Crespo, A., k Reach, 
W. T. 2006, arXiv:astro-ph/0604238 

Forster Schreiber, N. M., Roussel, H., Sauvage, M., k Charmandaris, V. 2004, A&A, 419, 501 

Furlan, E., Hartmann, L., Calvet, N., D'Alessio, P., Franco-Hernandez, R., Forrest, W. J., Watson, 
D. M., Uchida, K. L, Sargent, B., Green, J. D., Keller, L. D., k Herter, T. L. 2006, ApJS, 
165, 568 

Genzel, R., Lutz, D., Sturm, E., Egami, E., Kunze, D., Moorwood, A. F. M., Rigopoulou, D., 
Spoon, H. W. W., Sternberg, A., Tacconi-Garman, L. E., Tacconi, L., k Thatte, N. 1998, 
ApJ, 498, 579 

Helou, G., Roussel, H., Appleton, P., Frayer, D., Stolovy, S., Storrie-Lombardi, L., Hurt, R., 
Lowrance, P., Makovoz, D., Masci, F., Surace, J., Gordon, K. D., Alonso-Herrero, A., En- 
gelbracht, C. W., Misselt, K., Rieke, G., Rieke, M., Willner, S. P., Pahre, M., Ashby, 
M. L. N., Fazio, G. G., k Smith, H. A. 2004, ApJS, 154, 253 

Higdon, S. J., Higdon, J. L., k Marshall, J. 2006, ApJ, 640, 768 

Hogg, D. W., Tremonti, C. A., Blanton, M. R., Finkbeiner, D. P., Padmanabhan, N., Quintero, 
A. D., Schlegel, D. J., k Wherry, N. 2005, ApJ, 624, 162 

Houck, J. R., Charmandaris, V., Brandl, B. R., Weedman, D., Herter, T., Armus, L., Soifer, B. T., 
Bernard-Salas, J., Spoon, H. W. W., Devost, D., k Uchida, K. I. 2004, ApJS, 154, 211 



53 



Hudgins, D. M. & Allamandola, L. J. 1999, ApJ, 513, L69 

Hudgins, D. M. & Allamandola, L. J. 2005, in IAU Symposium, ed. D. C. Lis, G. A. Blake, & 
E. Herbst, Vol. 231, 443-454 

Hunt, L., Bianchi, S., & Maiolino, R. 2005, A&A, 434, 849 

Irwin, J. A. & Madden, S. C. 2006, A&A, 445, 123 

Joblin, C, Abergel, A., Bregman, J., D'Hendecourt, L., Heras, A. M., Jourdain de Muizon, M., 
Pech, C, & Tielens, A. G. G. M. 2000, in ESA SP-456: ISO Beyond the Peaks: The 2nd ISO 
Workshop on Analytical Spectroscopy, ed. A. Salama, M. F. Kessler, K. Leech, & B. Schulz, 
49-54 

Jones, A. P. & d'Hendecourt, L. 2000, A&A, 355, 1191 

Jones, A. P., Duley, W. W., & Williams, D. A. 1990, QJRAS, 31, 567 

Jura, M. 1987, in NATO ASIC Proc. 191: Polycyclic Aromatic Hydrocarbons and Astrophysics, 
ed. A. Leger, L. D'Hendecourt, & N. Boccara, 3-12 

Jura, M. 2003, ApJ, 582, 1032 

Jura, M., Bohac, C. J., Sargent, B., Forrest, W. J., Green, J., Watson, D. M., Sloan, G. C, 
Markwick-Kemper, F., Chen, C. H., & Najita, J. 2006, ApJ, 637, L45 

Kahanpaa, J., Mattila, K., Lehtinen, K., Leinert, C, & Lemke, D. 2003, A&A, 405, 999 

Kaneda, H., Onaka, T., & Sakon, I. 2005, ApJ, 632, L83 

Kawada, M., Shibai, H., Kaneda, H., & Nakagawa, T. 2004, in Optical, Infrared, and Millimeter 
Space Telescopes. Edited by Mather, John C. Proceedings of the SPIE, Volume 5487, pp. 
359-368 (2004)., ed. J. C. Mather, 359-368 

Kim, S.-H., Martin, P. G., & Hendry, P. D. 1994, ApJ, 422, 164 

Lagache, G., Dole, H., Puget, J.-L., Perez- Gonzalez, P. G., Le Floc'h, E., Rieke, G. H., Papovich, 
C, Egami, E., Alonso-Herrero, A., Engelbracht, C. W., Gordon, K. D., Misselt, K. A., & 
Morrison, J. E. 2004, ApJS, 154, 112 

Latter, W. B. 1991, ApJ, 377, 187 

Leger, A. & Puget, J. L. 1984, A&A, 137, L5 

Li, A. 2005, ApJ, 622, 965 

Li, A. & Draine, B. T. 2001a, ApJ, 554, 778 

-. 2001b, ApJ, 550, L213 



54 



-. 2002a, ApJ, 564, 803 
-. 2002b, ApJ, 572, 232 
-. 2002c, ApJ, 576, 762 

Li, A. & Greenberg, J. M. 1997, A&A, 323, 566 
-. 2002, ApJ, 577, 789 

Linsky, J. L., Draine, B. T., Moos, H. W., Jenkins, E. B., Wood, B. E., Oliveira, C, Blair, W. P., 
Friedman, S. D., Gry, C., Knauth, D., Lacour, S., Lehner, N., Redfield, S., Shull, J. M., 
Sonneborn, G., & Williger, G. 2006, ApJ, 647, 1106 

Lu, N., Helou, G., Werner, M. W., Dinerstein, H. L., Dale, D. A., Silbermann, N. A., Malhotra, S., 
Beichman, C. A., & Jarrett, T. H. 2003, ApJ, 588, 199 

Lutz, D., Valiante, E., Sturm, E., Genzel, R., Tacconi, L. J., Lehnert, M. D., Sternberg, A., & 
Baker, A. J. 2005, ApJ, 625, L83 

Madden, S. C., Galliano, F., Jones, A. P., & Sauvage, M. 2006, A&A, 446, 877 

Malloci, G., Joblin, C., & Mulas, G. 2006, in preparation, 000, 000 

Malloci, G., Mulas, G., & Joblin, C. 2004, A&A, 426, 105 

Mason, R. E., Wright, G. S., Adamson, A., & Pendleton, Y. 2006, astro-ph/0611374 
Mathis, J. S. 1996, ApJ, 472, 643 

Mathis, J. S., Mezger, P. G., & Panagia, N. 1983, A&A, 128, 212 
Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 
Mathis, J. S. & Whiffen, G. 1989, ApJ, 341, 808 

Mattioda, A. L., Allamandola, L. J., & Hudgins, D. M. 2005a, ApJ, 629, 1183 
Mattioda, A. L., Hudgins, D. M., & Allamandola, L. J. 2005b, ApJ, 629, 1188 
Mennella, V., Brucato, J. R., Colangeli, L., & Palumbo, P. 1999, ApJ, 524, L71 
Moutou, C., Leger, A., & D'Hendecourt, L. 1996, A&A, 310, 297 
Moutou, C., Sellgren, K., Verstraete, L., & Leger, A. 1999, A&A, 347, 949 
Mulas, G., Malloci, G., Joblin, C., & Toublanc, D. 2006, astro-ph/0606264, 000, 000 



55 



Onaka, T., Fujishiro, N., Ihara, O, Ishihara, D., Ita, Y., Kataza, H., Kim, W., Matsuhara, H., 
Matsumoto, T., Murakami, H., Sakon, I., Uemizu, K., Ueno, M., Wada, T., & Watarai, 
H. 2004, in Optical, Infrared, and Millimeter Space Telescopes. Edited by Mather, John C. 
Proceedings of the SPIE, Volume 5487, pp. 338-349 (2004)., ed. J. C. Mather, 338-349 

Onaka, T., Yamamura, I., Tanabe, T., Roellig, T. L., & Yuen, L. 1996, PASJ, 48, L59 

Pagani, L., Lequeux, J., Cesarsky, D., Donas, J., Milliard, B., Loinard, L., & Sauvage, M. 1999, 
A&A, 351, 447 

Papoular, R., Ellis, K., Guillois, O., Reynaud, O, & Nenner, I. 1993, J. Chem. Soc. Faraday Trans, 
89, 2289 

Peeters, E., Allamandola, L. J., Bauschlicher, Jr., C. W., Hudgins, D. M., Sandford, S. A., & 
Tielens, A. G. G. M. 2004, ApJ, 604, 252 

Peeters, E., Hony, S., Van Kerckhoven, O, Tielens, A. G. G. M., Allamandola, L. J., Hudgins, 
D. M., & Bauschlicher, C. W. 2002, A&A, 390, 1089 

Peeters, E., Tielens, A. G. G. M., Boogert, A. C. A., Hayward, T. L., & Allamandola, L. J. 2005, 
ApJ, 620, 774 

Piovan, L., Tantalo, R., & Chiosi, C. 2006, MNRAS, 366, 923 
Plante, S. & Sauvage, M. 2002, AJ, 124, 1995 

Popescu, C. O, Misiriotis, A., Kylafis, N. D., Tuffs, R. J., & Fischera, J. 2000, A&A, 362, 138 

Reach, W. T., Boulanger, F., Contursi, A., & Lequeux, J. 2000, A&A, 361, 895 

Roche, P. F., Aitken, D. K., Smith, C. H., & Ward, M. J. 1991, MNRAS, 248, 606 

Roelfsema, P. R., Cox, P., Tielens, A. G. G. M., Allamandola, L. J., Baluteau, J. -P., Barlow, M. J., 
Beintema, D., Boxhoorn, D. R., Cassinelli, J. P., Caux, E., Churchwell, E., Clegg, P. E., 
de Graauw, T., Heras, A. M., Huygen, R., van der Hucht, K. A., Hudgins, D. M., Kessler, 
M. F., Lim, T., & Sandford, S. A. 1996, A&A, 315, L289 

Rosenberg, J. L., Ashby, M. L. N., Salzer, J. J., & Huang, J.-S. 2006, ApJ, 636, 742 

Sakata, A., Wada, S., Onaka, T., & Tokunaga, A. T. 1987, ApJ, 320, L63 

Sakon, I., Onaka, T., Ishihara, D., Ootsubo, T., Yamamura, I., Tanabe, T., & Roellig, T. L. 2004, 
ApJ, 609, 203 

Salama, F., Bakes, E. L. O., Allamandola, L. J., & Tielens, A. G. G. M. 1996, ApJ, 458, 621 
Sellgren, K., Luan, L., & Werner, M. W. 1990, ApJ, 359, 384 



56 



Siebenmorgen, R. & Kriigel, E. 1992, A&A, 259, 614 

Silva, L., Granato, G. L., Bressan, A., & Danese, L. 1998, ApJ, 509 

Smith, J. D. T., Dale, D. A., Armus, L., Draine, B. T., Hollenbach, D. J., Roussel, H., Helou, G., 
Kennicutt, R. C, Li, A., Bendo, G. J., Calzetti, D., Engelbracht, C. W., Gordon, K. D., 
Jarrett, T. H., Kewley, L., Leitherer, C, Malhotra, S., Meyer, M. J., Murphy, E. J., Regan, 
M. W., Rieke, G. H., Rieke, M. J., Thornley, M. D., Walter, F., & Wolfire, M. G. 2004a, 
ApJS, 154, 199 

Smith, J. D. T., Draine, B. T., Dale, D. A., Moustakas, J., Kennicutt, R. C, Helou, G., Armus, 
L., Roussel, H., Sheth, K., Bendo, G. J., Buckalew, B., Calzetti, D., Engelbracht, C. W., 
Gordon, K. D., Hollenbach, D. J., Li, A., Malhotra, S., Murphy, E. J., & Walter, F. 2006, 
ApJ, submitted, 000, 000 

Smith, T. L., Clayton, G. C, & Valencic, L. 2004b, AJ, 128, 357 

Snow, T. P. 2000, J. Geophys. Res., 105, 10239 

Snow, T. P. & Witt, A. N. 1996, ApJ, 468, L65 

Speck, A. K. & Barlow, M. J. 1997, Ap&SS, 251, 115 

Tanabe, T., Nakada, Y., Kamijo, F., & Sakata, A. 1983, PASJ, 35, 397 

Tanaka, M., Matsumoto, T., Murakami, H., Kawada, M., Noda, M., & Matsuura, S. 1996, PASJ, 
48, L53 

Thuan, T. X., Sauvage, M., & Madden, S. 1999, ApJ, 516, 783 

Tielens, A. G. G. M. 2005, The Physics and Chemistry of the Interstellar Medium (The Physics and 
Chemistry of the Interstellar Medium, by A. G. G. M. Tielens, pp. . ISBN 0521826349. Cam- 
bridge, UK: Cambridge University Press, 2005.) 

Tielens, A. G. G. M., Hony, S., van Kerckhoven, C, & Peeters, E. 1999, in ESA SP-427: The 
Universe as Seen by ISO, ed. P. Cox & M. Kessler, 579-587 

Tuffs, R. J., Popescu, C. C, Volk, H. J., Kylafis, N. D., & Dopita, M. A. 2004, A&A, 419, 821 

Uchida, K. I., Sellgren, K., & Werner, M. 1998, ApJ, 493, L109 

Uchida, K. I., Sellgren, K., Werner, M. W., & Houdashelt, M. L. 2000, ApJ, 530, 817 

van Diedenhoven, B., Peeters, E., Van Kerckhoven, C, Hony, S., Hudgins, D. M.and Allamandola, 
L. J., & Tielens, A. G. G. M. 2004, ApJ, 611, 928 

Verstraete, L., Pech, C, Moutou, C, Sellgren, K., Wright, C. M., Giard, M., Leger, A., Timmer- 
mann, R., & Drapatz, S. 2001, A&A, 372, 981 



57 



Verstraete, L., Puget, J. L., Falgarone, E., Drapatz, S., Wright, C. M., k Timmermann, R. 1996, 
A&A, 315, L337 

Webster, A. 1993, MNRAS, 264, 121 

Weingartner, J. C. k Draine, B. T. 2001a, ApJ, 548, 296 

-. 2001b, ApJS, 134, 263 

Werner, M. W., Roellig, T. L., Low, F. J., Rieke, G. H., Rieke, M., Hoffmann, W. F., Young, E., 
Houck, J. R., Brandl, B., Fazio, G. G., Hora, J. L., Gehrz, R. D., Helou, G., Soifer, B. T., 
Stauffer, J., Keene, J., Eisenhardt, P., Gallagher, D., Gautier, T. N., Irace, W., Lawrence, 
C. R., Simmons, L., Van Cleve, J. E., Jura, M., Wright, E. L., k Cruikshank, D. P. 2004a, 
ApJS, 154, 1 

Werner, M. W., Uchida, K. I., Sellgren, K., Marengo, M., Gordon, K. D., Morris, P. W., Houck, 
J. R., k Stansberry, J. A. 2004b, ApJS, 154, 309 

Witt, A. N. k Gordon, K. D. 1996, ApJ, 463, 681 

Wu, Y., Charmandaris, V., Hao, L., Brandl, B. R., Bernard- S alas, J., Spoon, H. W. W., k Houck, 
J. R. 2006, ApJ, 639, 157 

Yan, L., Chary, R., Armus, L., Teplitz, H., Helou, G., Frayer, D., Fadda, D., Surace, J., k Choi, 
P. 2005, ApJ, 628, 604 

Zhou, Z., Sfeir, M. Y., Zhang, L., Hybertsen, M. S., Steigerwald, M., k Brus, L. 2006, ApJ, 638, 
L105 

Zubko, V., Dwek, E., k Arendt, R. G. 2004, ApJS, 152, 211 



This preprint was prepared with the AAS IATjrjX macros v5.2. 



